Empirical Relationships for 
Peak Ground Acceleration in Indian Earthquakes 




A Thesis Submitted 

In Partial Fulfillment of the Requirements 
for the Degree of 

MASTER OF TECHNOLOGY 


by 

Roshan A. D. 
Roll number 9710328 



DEPARTMENT OF CIVIL ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY KANPUR 

Aueust 2000 




jwl'rfT 

sTrsYfi-i^V 

3r5rif% spo A-. m2«5 — 
TH’ 




CERTIFICATE 


It is certified that the work contained in this 
Relationships for Peak Ground Acceleration in Indian Earthquakes” by, Mr. Roshan 
A.D. (Roll Number: 9710328), has been carried out under our supervision and this work 
has not been submitted elsewhere for a degree. 



entitled “Empirical 


(Dr. Pr|i^rb. Basu) 

Head, Civil & Structural Engineering Division 
Atomic Energy Regulatory Board 
Niyamak Bhavan, Mumbai 


(Dr. SudhifE"jain) 
Professor 

Department of Civil Engineering 
Indian Institute of Technology, Kanpur 


ABSTRACT 


Himalayan region in India is one of the most seismically active regions in the world. 
About 120 recordings from Strong Motion Accelerographs (SMAs) are now available from 
the recent Himalayan earthquakes in the magnitude {M) range of 5.2 to 7.3 with the epicentral 
distances {R) ranging from 3 to 350kms. About 260 records from Structural Response 
Recorders (SRRs) from 5 Himalayan earthquakes are also available (M 5.5-7.3; R: 4- 
770kms). SRRs record spectral acceleration at natural periods of 0.4sec, 0.75sec and 1.25sec, 
for damping ratios of 5% and 10% of critical. The strong motion database available for Indian 
earthquakes can be enhanced by making use of records from SRRs. 

The compatibility of responses from SRRs and SMAs from Indian earthquakes has 
been studied and the peak ground acceleration has been estimated from SRR recordings. The 
estimates of peak ground acceleration from SRRs are in line with those from SMAs. The 
peak ground accelerations in Indian earthquakes are compared with those predicted by the 
existing attenuation relationships. It is seen that the relationships provide reasonable 
estimates for the earthquakes in Central Himalayas, but underestimate the motion for 
earthquakes in North-East India and Indo-Gangetic plains. It is seen that based on attenuation 
characteristics, the Himalayan earthquakes can be grouped into four categories; (a) events of 
Central Himalayas, (b) non-subduction earthquakes of North-East India, (c) subduction 
events of North-East India, and (d) Earthquakes in Indo-Gangetic plains. 

Separate attenuation relationships for peak ground acceleration are derived for 
earthquakes in Central Himalayas, Indo-Gangetic plains. North East Indian Non-subduction 
zone earthquakes and North East Indian subduction zone earthquakes. 
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Chapter 1 

Introduction 

1.0 General 

The Probabilistic Risk Assessment (PRA) of a Nuclear Power Plant (NPP) is 
conducted to develop an appreciation of plant behavior, understand most likely accident 
sequences, identify the dominant risk contributors etc. The contribution to the total plant 
risk can be divided into two parts; (1) Risk due to internal events, (2) Risk due to 
external events. The risk due to seismic activity is considered to be a major contributor in 
the total risk due to external events. This is because of the unique ability of an 
earthquake to initiate an accident and simultaneously fail a number of otherwise 
redundant components required for mitigating the accident. The elements of a seismic 
risk analysis can be identified as analyses of: 

1 . the seismic hazard at site 

2. response of plant systems and structures 

3. component fragilities 

4. accident sequences 

5. consequences 

1.1 Seismic hazard analysis 

The seismic hazard assessment can be carried out either deterministically or 
probabilistically. 

1.1.1 Deterministic Seismic Hazard Assessment 

The most straightforward approach to seismic Hazard assessment is 
deterministic. A Deterministic Seismic Hazard Assessment (DSHA) is based on 
identifying scenario events and determining the resulting ground motions at site. The 
basic steps of DSHA are illustrated in Figure 1.1. The size of the design earthquake may 



Chapter 1:: Introduction 2 

be determined from the length of faults, assuming that a certain portion will rupture and 


estimating the magnitude from empirical relationships. When this is not possible, it is 
common practice to estimate the maximum credible earthquake magnitude, not from the 
physical dimensions of the fault but by simply adding an increment to the largest 
earthquake in the seismic catalogue [Bommer, 1999], 

Having fixed the design magmtude for each source and corresponding distance 
from source, it is a simple matter to calculate the resulting ground motion at site with the 
help of attenuation relationships. The most severe case is then chosen as design 
earthquake. But a major difficulty with DSHA is that the level of conservatism of the 
design event is unknown. 

1.1.2 Probabilistic Seismic Hazard Assessment 

The problems associated with defining design earthquake are bypassed in Probabilistic 
Seismic Hazard Assessment (PSHA) (Figure 1.2). 

The major elements of the hazard analysis are: 

1. Identification of the sources of earthquakes such as faults and seismotectonic 
provinces. The assumption is made that earthquakes have an equal probability of 
occurring at any location within the source, so it is not necessary to fix the location of 
the design earthquake. 

2. Evaluation of the earthquake history of the region to assess the frequencies of 
occurrence of earthquakes of different magnitudes. 

3. Development of an attenuation relationship to estimate the intensity of earthquake 
induced ground motion (e.g., peak ground acceleration) at the site. 

4. Integration of all these information to generate the frequencies with which different 
values of selected ground motion parameter would exceed. 
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Thus, the attenuation of ground motion parameter depending on the distance from the 

source and the magnitude of the earthquake is one of the important parameters used in 
seismic hazard evaluation of site. 

The ground motion parameter predicted by the attenuation relation can be peak 
ground acceleration, peak gromid velocity, peak groimd displacement, intensity of 
shaking, or shape of the response spectrum. 

The attenuation relationships developed are constantly updated based on the 
availability of new data for the particular region. In Indian subcontinent, the attenuation 
relationships commonly used for major projects are based on the data from Western 
United States of America (e.g., Chandrasekaran and Das, 1992). In recent years, strong 
motion data has been recorded -by the strong motion instmmentation during several 
earthquakes in Northern India (e.g., Chandrasekaran and Das, 1990; Chandrasekaran and 
Das, 1992). 

1.2 Scope of Study 

A methodology is developed to calculate the peak groxmd acceleration from the 
response of stmctural response recorders. Based on this, a new attenuation relationship 
for peak ground acceleration is proposed using the Strong Motion Accelerograms and 
Structural Response Recordings of earthquakes experienced in the Himalayan region 
from 1986 to 1999. 

1.3 Organisation 

The work has been presented in five chapters including this one. In Chapter 2, 
various parameters, which control the attenuation of acceleration are discussed along 
with different solution schemes available. A review of various attenuation models 
available in the literature for different regions of the world is also presented. Chapter 3 
discusses the various types of strong motion records available for Indian earthquakes and 
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formulates methods of extraction of peak ground acceleration from these records. A brief 

comparison of the observed peak ground acceleration with some of the published 

attenuation equations is presented in Chapter 4. A set of new attenuation relationships is 

also developed for the Himalayan region. A brief summary and conclusions of ±e study 

are presented in Chapter 5. The database used in the present study is given in Appendices 

A and B. 

• • • 
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Chapter 2 

Review of Existing Attenuation Relationships 

2.0 Introduction 

Seismic hazard assessment may be performed deterministically and 
probabilistically. Both approaches require the use of ground motion attenuation models 
for prediction of ground motion parameters. Most of the attenuation models are based on 
statistical analysis of recorded ground motions, which are updated as new strong motion 
data becomes available. This chapter reviews some of the attenuation relationships 
available in the literature for calculation of peak ground acceleration and spectral 
ordinates of response spectra. 

2.1 Estimation of Ground Motion 

The ground motion at a site depends on magnitude of the earthquake (M), distance 
{R\ site conditions (5), type of faulting {F) etc. Different source-to-site distance 
measurements used by different researchers are given in Figure 2.1 and 2.2. Different 
type of site classification schemes ranging from qualitative descriptions of the base 
surface material [Calliot and Bard, 1993] to very quantitative definitions based on shear 
wave velocity [Boore et al., 1997] are used. 

With the availability of more strong motion data, separate relationships have been 
developed for different categories of regional ground motion characteristics. These 
include shallow crustal earthquakes in active tectonic regimes, shallow crustal 
earthquakes in stable continental regions, subduction zone earthquakes, and earthquakes 
in extensional regimes. Regardless of the tectonic regime, rupture directivity also affects 
ground motion attenuation relationships [Somerville et al, 1997] 

Attenuation relationships are available for predicting not only the peak ground 



Chapter 2:: Review of Existing Attenuation . . . 


8 

acceleration, but also the spectral ordinates of the response spectrum. In all the available 
relationships, the ordinates of response spectrum are calculated for 5% of critical 
damping. 

2.2 Attenuation Relationships for Peak Ground Acceleration (PGA\ 

2.2.1 Nature of PGA Values 

The strong motion accelerographs provide three components of ground motion: 
two horizontal (longitudinal and transverse) and one vertical. For earthquake resistant 
design, one is usually more interested in estimation of horizontal acceleration. There are a 
number of methods for defining even a simple parameter like PGA since each 
accelerogram has two horizontal components.: 

1. The two horizontal components are considered as separate and independent records 
during the regression analysis [Crouse, 1990]. This method corresponds to a 
randomly oriented component and hence allows statistical treatment of the variation 
in the orientation of the shaking with respect to the axes of the structure. But, this 
procedure may lead to erroneous results due to high correlation that exists between 
these two components [Campbell, 1985]. 

2. The largest of the two horizontal components is considered for the regression analysis 
[Joyner and Boore, 1981; Anders etal., 1990; Tento et al., 1992; Ambraseys, 1995; 
Loh et al., 2000]. It has been reported that on an average, the larger value of the peak 
acceleration exceeds the mean value of the acceleration by about 13 percent [Joyner 
and Boore, 1981]. 

3. The arithmetic mean of the two horizontal components is considered in the regression 
analysis [Fukushima and Tanaka, 1990; Boore et al, 1997; Campbell, 1997; Sadigh 
et al, 1997]. The regression analysis using this method results in almost identical 
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predicted values as those obtained by using both components of PGA, but artificially 


reduces scatter in the regression [Bommer, 1999]. 

4. The geometric mean of the two components is used in the regression analysis [Young 
etal, 1997]. 

2.2.2 Propagation Parameters 

The ground motion attenuation is dependent on the propagation parameters of the 
medium as well as the source parameters. Propagation parameters characterize the wave 
scattering, geometrical attenuation, and inelastic attenuation of the ground motion as it 
travels from the source to the site. 

2.Z2.1 Distance 

The variable that is generally used to characterize the propagation parameter is 
distance. The earthquake rupture can extend over tens to hundreds of kilometers, and a 
number of distance measures have come into use. The various distance measures used in 
the literature are (Figure 2.1): (a) epicentral distance (i2e), (b) hypocentral distance {Rkyp^, 
(c) distance to zone of energy release (i?z), (d) closest distance to rupture (e) 

closest distance to surface projection of rupture {rj^. These are briefly described below. 
Epicentral distance (R^: 

It is the distance on the ground surface between an observer or site and the 
epicenter. Many relationships derived prior to 1985 were based on epicentral distance. 
Ambraseys (1995) has also used the epicentral distance because only a few earthquakes 
in the database had the information about associated surface faulting. When the 
information about the fault dimensions were available, the shortest distance to surface 
projection of the fault {rj^ was used as the distance measure and in the case of events 
where the fault parameters were not known, rjb was assumed equal to epicentral distance. 
Hypocentral distance (Rhypo)- 

Atkinson and Boore (1997a) have used hypocentral distance as the distance 
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parameter for modeling the ground motion attenuation in the Cascadia region. The study 
conducted by Joyner and Boore (1981) suggest that a fault distance based measure is 
more appropriate in arriving at an attenuation relationship compared to the one based on 
distance to focus or epicenter. It was observed in their study that the closest distance to 
rupture zone is most suitable representation of the distance especially in the near source 
region. However, due to uncertainties in the depth estimation, this factor can be replaced 
by distance to surface projection of fault [Calliot and Bard, 1993]. Their study compared 
the results of regression using hypocentral distance with those using distance to surface 
projection of fault as the distance measure. It was foimd that the different distance 
definitions show very little variation in the statistical parameters of the results especially 
for events with magnitude less than 6.5. 

Distance to zone of energy release (R^\ 

During the fault rupture, the maximum amount of energy released may be 
originating from a small zone in the fault plane. The distance to this zone can be used as a 
distance measure in the attenuation relationship. During the development of attenuation 
relationship for subduction zone earthquakes, Crouse (1991) used the distance to center 
of energy release as the distance measure. Tento et al., (1992) considered Rz to be 

equivalent to i? = -yjd^ + Where, d is the closest distance to fault and h is hypocentral 
depth. In the study, the use of i? as the distance measure resulted in lesser scatter of the 
data compared to that using J as the distance measure. 

Closest distance to plane of rupture (Rrup)- 

The shortest distance between the site and fault rupture was used as the distance 
measure in the development of some attenuation relationships [Fukushima and Tanaka, 
1990; Abrahamson and Silva, 1997; Sadigh et al, 1997]. Based on the recordings 
obtained from Chi-Chi earthquake of Taiwan, 21 September 1999, Loh et al., (2000) 
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developed an attenuation relationship with the shortest distance to plane of rupture being 

us^ed as the distance measure. Campbell (1994) used distance to seismogenic rupture 
{Rseis) tti the derivation of attenuation relationship. The seismogenic rupture zone was 
determined from the location of the surface fault rupture, the spatial distribution of 
aftershocks, earthquake modeling studies, regional crustal velocity profiles, and geologic 
data. This distance measure was used under the assumption that the fault mpture wi thin 
the upper sediments and within 2 to 4 km of the fault is primarily non-seismogenic and 
does not contribute significantly to the recorded ground motion at frequencies of 
engineering interest. 

Closest distance to surface projection of fault rupture (rjb)‘. 

Joyner and Boore (1981) have used closest distance to surface projection of fault 
rupture as the distance measure. Therefore, this distance is also known as Joyner-Boore 
distance (rjb). The epicentral distance or distance to center of rupture was not used in their 
analysis because during the Parkfield (1966) and Imperial Valley (1979) earthquakes the 
recording sites were located close to rupture but far from epicenter, and rupture center 
recorded high values of acceleration. Though the distance to point on the rapture would 
have been a more appropriate measure of distance, it was not used in their relationship 
because of the difficulties involved in determination of that point. Many of the recent 
relationships [Spudich et al, 1997; Toro et al., 1997] also use this as the distance 
measure in their calculations. 

For the sites located several source dimensions from the earthquake (far source 
records), there is little difference between the results obtained using different distance 
measures. But for shorter distances (near source records), the difference in results using 
different distance measures become more significant. In the near source region, where 
predictions are of greatest concern, the use of epicentral or hypocentral distance leads to 
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considerably greater scatter in estimates of strong ground motion. Use of a fault-distance 

.measure can lead to biased predictions, especially if the strong-motion stations included 
in the analysis have a non-random distribution around the fault, or if the strong motions 
come from a localized source (or sources) then on the fault [Campbell, 1985]. The 
concept of a distance to localized source is represented by M3 in Figure 2.2. 

If the strong motion radiates from small areas of the fault rupture surface, then a 
fault-distance measure would tend to underestimate the actual distance to these localized 
sources. It may be possible to identify these areas for some past earthquakes, but it is 
virtually impossible to anticipate their locations during future events. Because of this, 
most earthquake scenarios, whether for probabilistic or deterministic applications, use the 
closest approach to the fault, tectonic structure, or earthquake rupture as the 
representative distance of the hypothesized earthquake. If an analysis hypothesizes 
earthquake sources to be equally distributed along a fault or within an area, then 
epicentral distance, hypocentral distance, or distance to the energy center would be the 
more appropriate measure to use. In this case, attenuation relations in terms of fault 
distance will underestimate the true ground motions. This is because, for a fault under 
consideration the distance measures like epicentral distance, hypocentral distance or 
distance to energy center will vary based on the originating point of rupture whereas 
closest distance to surface projection of fault will remain same. 

2.2.2.2 Site Geology 

Site parameters in a model are related to geologic descriptions of recording 
stations. The recorded ground motion may be affected by the nature of soil on which the 
recording station is situated, A large amplification (as much as a factor of two) in 
accelerations was observed to be associated with shallow deposits for sites located near 
the source of small to moderate earthquakes [Campbell, 1985]. The classification of these 



Chapter 2:: Review of Existing Attenuation , . . 13 

sites as rock can significantly increase estimates of strong ground motion if enough of 

these sites are included in the analysis. The different classification schemes used by 

researchers divide the sites into hard rock, soft rock, shallow alluvium, deep alluvium. 

stiff soils, deep cohesionless soils, soft soils, etc. The classification is based on the depth 

of the soil deposit, age of the deposit, shear wave velocity of the medium, etc. 

The soil response can be modeled either as a function of peak ground acceleration 
or by using an amplification factor, which predicts the soil site response based on the 
corresponding motion in the rock. Thus, a single attenuation relationship can be used for 
both soil and rocky strata. It is also possible to use separate fimctional fo rms for soil and 
rock resulting in separate regression analysis for soil and rock site recordings. Yoimgs et 
al. (1997) classified the sites into rock, shallow soil, and deq) soil. Separate relationships 
were provided for rocky and deep soil sites. For an earthquake of magnitude 5.0, the 
predicted accelerations on deep soil sites were 50% higher than that for rocky sites. This 
difference increased to 65% for an earthquake of magnitude 7.0. Abrahamson and Silva 
(1997) classified the sites into rock and deep soil, which was based upon geomatrix soil 
classification. The amplification factor for soil was predicted as a nonlinear function of 
PGA in rocky strata. The predicted relationship showed an higher PGA in rocky sites 
compared to soil sites up to a distance of 40 km. After 40km, predicted PGA was higher 
for soil sites. The difference in predicted PGA for soil sites varied from -30% at 5km to 
+30% at 300km with respect to rocky sites for an earthquake of magnitude 7.0 

Boore et al. (1997) used shear wave velocity of upper 30m of the crust to 
characterize the sites into 5 groups namely NEHRP site class B, class C, class D, rock 
and soil. When compared with the response predicted by hJEHRP class B; NEHRP site 
class C, class D, rock and soil predicted the PGAs, which were higher by about 30%, 
71%, 22% and 58% respectively. Campbell (1997) classified the sites into alluvium or 
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firm soil with depth of layer greater than 10m, soft rock, and hard rock. As per this 
classification, for an earthquake of magnitude 7.0, acceleration predicted for soil site was 
96% higher than that of the rocky sites. This was reduced to 35% in the case of sites with 
soft rocky strata. The site dependent amplification factor used by Campbell (1997) was a 
fimction of both magnitude and distance. 

2.2.3 Effect of Magnitude 

The parameter mostly used to characterize the earthquake size in strong motion 
attenuation relationships is earthquake magnitude. It is the only source parameter 
routinely reported by seismographic networks. Earthquake magnitude although routinely 
reported and universally used as a measure of earthquake size has a few limitations. The 
variety of magnitude scales that exist can lead to confusion in comparing various 
predictions. There is also a clear tendency for all scales, except moment magnitude, to 
saturate as the size of the earthquake increases (Figure 2.3). 

Currently, the following different types of magnitude scales are being used 

1. Richter Local Magnitude {Mi): This magnitude scale is based on the maximum trace 
amplitude of the Wood- Anderson seismometer. Commonly used for shallow, local 
earthquakes. 

2. Surface Wave magnitude {Ms): The amplitude of Rayleigh waves is used to calculate 
surface wave magnitude. The surface wave magnitude is commonly used to describe 
the size of shallow focus, distant, moderate to large earthquakes. 

3. Body wave magnitude ( ): For deep focus earthquakes, surface waves are often too 

small to permit reliable evaluation of the surface wave magnitude. The body wave 
magnitude scale is based on the amplitude of the first few cycles of p-waves (usually 
with a period of Isec) which are not strongly influenced by the focal depth. 



Chapter 2:: Review of Existing Attenuation . . . 15 

4. Moment magnitude {Mw)'- The ground shaking characteristics does not increase at the 

same rate at which the total energy released increases during an earthquake. For 

strong earthquakes, the measured ground shaking characteristics become less 

sensitive to the size of the earthquake, than for smaller earthquakes. The body wave 

and Richter local magnitudes saturate at magnitudes of 6.0 to 7.0 and the surface 

wave magnitude saturates at about Ms = 8.0. In order to overcome this drawback, a 

magnitude scale based on the seismic moment was developed. This magnitude scale 

is known as moment magnitude. Figure 2.3 gives the comparison of different 

magnitude scales being used worldwide with moment magnitude. 

An important element in the choice of scale is related to the specification of the 
magnitude of a future earthquake. Generally, surface wave magnitudes are not reliably 
determined for magnitudes less than 6.0 and because of saturation become relatively 

independent of earthquake size for magnitudes near 7.0. 

Most of the recent attenuation relationships [Boore et al, 1991 \ Youngs et al, 
1997; Joyner and Boore, 1981; Crouse, 1991] make use of moment magnitude. Since Mi 
is approximately equal to Mw for Mi<6.5 and Ms is approximately equal to Mw for Ms 
ranging from 6.0 to 8.0, Campbell (1997) made use of Mi and Ms magnitudes in the 
interval as specified above considering them as equal to Mw- Ambraseys (1995) while 
developing attenuation relationships for European earthquakes, used the events which 
had their surface wave magnitude (Ms) between 4.0 and 7.3. Calliot (1993) made use of 
Ms if (Ms and Mi) > 6.0 and Mi if (Ms or Mi) < 6.0 as the magnitude term. Fukushima 
and Tanaka (1990) used Mj, the magnitude defined by Japan Meteorological Agency for 
earthquakes with My >6.0. 
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2.2.4 Source Mechanism 

Different tectonic environments give rise to different ground motion attenuation 
characteristics. Recent studies [Campbell, 1997] on ground motions show that normal- 
faulting earthquakes located in extensional stress regimes are associated with lower 
ground motions than either strike-slip or reverse faulting earthquakes located in 
compressional stress regimes. For distances greater than SOkm from the earthquake 
rupture, ground motions from subduction zone earthquakes are substantially larger than 
those from shallow crustai earthquakes [Crouse, 1991]. 

The distinction between ground motions from strike slip and reverse faults has 
become common in recent relationships [Campbell and Bozorgnia 1994; Abrahamson 
and Silva, 1997]. The amplification due to style of faulting can be modeled as a value 
independent of both source-to-site distance and magnitude [Boore et al, 1997], as a 
function of only magnitude [Abrahamson and Silva, 1997] or as a function of the distance 
measure and magnitude [Campbell, 1997]. Boore et al. (1997) predicted an increase in 
acceleration by 21% in the case of reverse faults compared to strike-slip faults. As per 
Abrahamson and Silva (1997), the increase was 30% for reverse faults and 15% for 
oblique faults compared to strike slip events with a magnitude of 7.0. 

2.2.5 Method of Regression Analysis 

Different methods of regression analysis are being used by researchers to arrive at 
the attenuation relationships: 

2.2. 5.1 Method of Least Squares 

In the past, most of the regression analyses for curve fitting were performed using 
a least-squares procedure. This procedure minimizes the sum of the squared errors 
n 

, where y is the predicted value of y,- and w,- is the weight assigned to each 

j=l 

observation y,-. Ambraseys (1995) used this method for the development of attenuation 
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relationship for European earthquake data. If the equation is linear with respect to the 

coefficients to be determined, then standard linear least-squares procedure can be used. 
Otherwise non-linear methods of analysis should be used. Campbell (1997) and Crouse 
(1991) used generalised non-linear least square analysis to derive the attenuation 
relationships. Stage-wise regression procedures are useful to decide the important 
variables controlling the resultant parameter. If the model is non-linear, distribution of the 
nonlinear coefficients must be developed empirically using Monte-Carlo simulation 
[Campbell, 1985]. 

Biased estimates of the coefficients will be obtained if many recordings from a 
few earthquakes or recording sites dominate the data. This bias can be reduced by 
restricting the data sample to no more than a certain number of recordings from a given 
earthquake for a given site, and with the use of weighted regression procedures to 
equalise the impact of the recordings from individual earthquakes or from specified 
ranges of magnitude and distance [Caillot and Bard, 1993]. In order to increase the 
impact of larger values of Y, some researchers have carried out the regression analysis on 
Y rather than ln7 [Bolt and Abrahamson, 1982]. 

2.2.5.2 Two Step Regression Analysis 

Two-step regression procedure de-couples the determination of magnitude 
dependence from determination of distance dependence. This procedure is also known as 
two-step stratified regression analysis. The variables considered in regression equations 
usually can take values over some continuous range. Occasionally a factor is to be 
introduced, which has two or more distinct levels. For example, data may arise from three 
machines, or two factories etc. In such case we cannot set up a continuous scale for the 
variable “machine”, or “factory”. We must assign these variables some levels in order to 
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take account of the feet that these variables may have separate and distinct effects on the 

response. Such variables are usually-called dummy variables. 

Two-step regression analysis introduces a dummy variable Ei for each earthquake. 

This variable, £,■ will take a value equal to ‘1’ for all the records obtained from the z* 

event and ‘0’ for the records from all other events. This procedure de-couples the 

determination of magnitude dependence from determination of distance dependence. In 

the first stage, the regression is carried out using the form of equation as given below. 

N 

log y = R + b 2 R + b^S (2.1) 

1=1 

where N is the total number of earthquakes in the data sample, Ei is the dummy variable 
whose value is equal to ‘1’ for earthquake ‘i’ and ‘0’ otherwise. The values a,-, bi, b 2 and 
bs are determined using method of least squares. After a,- values are determined, they are 
used to fit a first or second order polynomial by least squares to represent the magnitude 
dependence. 

ai=a+jM+]M^ (2.2) 

The advantages of this method are [Joyner and Boore, 1981]: 

1. The data from a single earthquake is typically recorded over a limited range of 
distance. If the regression analyses were done simultaneously including the 
magnitude and distance terms, errors in measuring magnitude would affect the 
distance coefficient obtained from regression and vice versa. 

2. This method causes each earthquake to have the same weight in determining the 
magnitude dependence and each recording to have the same weight in determining 
the distance dependence. 

The standard error exy of the residuals from the regression is given by 



where is the standard error from the first stage of regression and is the standard 

error obtained from the second stage of regression. Many of the recent attenuation 
relationships are based on two-stage methods e.g., Joyner and Boore (1981); Fukushima 
and Tanaka (1990); Anders et al. (1990), Tento et.al (1992); Galliot and Bard (1993); 
Boore et al. (1997). Fukushima and Tanaka (1990) showed that for the data set 
considered for regression analysis, one-stage ordinary least square results (non-stratified 
regression analysis) were seriously erroneous. The error was attributed to the strong 
correlation between magmtude and distance and the resulting trade-off between 
magnitude dependence and distance dependence. The correct distance dependence, given 
by the two-stage method and verified by analyzing individual earthquakes separately, 
showed a much stronger decay of PGA with distance than that predicted by the one-stage 
ordinary least squares method. The model used for the analysis was 
log(y) = aM - b logi? + c . The distance coefficients were calculated for individual events 
as well as by one stage and two stage methods, and it was shown that the two-stage 
regression was giving closer value to the average distance coefficient of the data (Figure 
2.4). 

Campbell (1985) and Caillot and Bard (1993) reported that the two-step 
regression procedure precludes the optimum fit of the data and gives only a local minima 
of the regression surface. Moreover, it fails to account for the magnitude saturation at 
near source distances. The analysis of the data obtained from European earthquakes 
showed that for the database under consideration, one-stage and two-stage methods 
produced almost same shape of attenuation curve (Figure 2.5(a), 2.5(b)) [Ambraseys, 
1995]. The model used to derive the particular attenuation relationship was 
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(2.4) 

2.2.5.3 Method of Maximum Likelihood 

The method of maximum likelihood is a general method of finding estimated 
(fitted) values of parameters. The underlying concept is that the best estimate of a 
parameter is that giving the highest probability that the observed set of measurements 
will be obtained. The estimation process involves considering the observed data values as 
constants and the parameter to be estimated as a variable, and then using differentiation 
to find the value of the parameter that maximizes the likelihood fiinction. 

Let / (x, G) be the density function of the random variable x and let (^be the parameter of 
the density function. If we observe a random sample of XiJC 2 ,..JCn, then the maximum 
likelihood estimator of ^is that value which maximises the density function f(x,G). The 
maximum likelihood estimator has several desirable properties 

a) The estimator is efficient in the sense that there is no estimator with smaller variance. 

b) The estimator approaches the true population parameter asymptotically as the number 
of observations increases. The maximum likelihood method works best for large 
samples, where it tends to produce estimators with the smallest possible variance. 

The maximum likelihood estimate for the slope and intercept in simple linear 
regression is the same as the least squares estimates when the underlying distribution for 
Y is normal [e.g., Gujarati, 1995]. In general, however, the maximum likelihood and least 
square estimates need not be the same. Method of maximum likelihood can be used to 
calculate the estimated value of a coefficient in either one-stage or two-stage regression 
analysis instead of using least square methods. Abrahamson and Silva (1997) and Boore 
et al. (1997) have used maximum likelihood method to develop the attenuation 
relationships. Joyner and Boore (1993) reported that using Fukushima and Tanaka’s data, 
the one-stage maximum likelihood method gave essentially same results as the two-stage 
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method. But, these results did not match with the parameters obtained using ordinary 

least square method (Table 2.1). 

2.2.6 Validity of Result and Treatment of Uncertainty 

2.2.6. 1 Analysis of Residuals 

Different methods for regression analyses are derived based on the assumption 
that the random error term Sf is normally distributed. This error term is represented by the 
residuals obtained after regression. The normality assumption of the residuals is checked 
by y? [Chi-Square] test. The adequacy of the model with respect to the fitted variables 
are then checked by plotting the residuals against the independent variables M, R, etc. If 
these plots do not show any apparent trend, the model can be considered adequate. 

2.2. 6.2 Multiple Correlation Coefficient 

The goodness of the fit can be evaluated by the multiple correlation coefficient. 
The expression for multiple correlation coefficient is [e.g., Caillot and Bard, 1993]: 



where w,- is the weight assigned to each observation, y,- is the predicted value of the 
parameter y is the mean of observations . Rc varies from 0 to 1 (perfect correlation). 
The proportion of the variability of y/ ’s explained by the regression is given by R^ . If 
R^ is less than 0, the regression form used is unable to explain the variation in the data. 
Caillot and Bard (1993) reported that the value of Rc varied from 0.4 to 0.5 when two- 
stage regression analysis was used to predict both PGA and spectral acceleration. 
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2.2.63 Uncertainty in Prediction 

J[n addition to median ground motion estimates, the uncertaintv (standard 
deviation) is also important for seismic hazard analyses. The uncertainty in the prediction 
of the strong motion parameters is represented with the help of standard deviation ( cr) of 
the predicted variables. The predicted dependence of standard deviation is different for 
different models, and also depends on the data set. The standard deviation is reported as a 
constant value or as a value dependant on the predicted amplitude or dependant on the 
magnitude of the earthquake predicted. Usually a linear relationship with respect to the 
magnitude or amplitude of motion is used to denote the uncertainty in the prediction of 
ground motion [Campbell, 1997; Sadigh et al., 1997]. The attenuation models having 
magnitude dependant standard deviation predicts smaller standard deviations for the 
events with larger magnitudes, than smaller magnitude events. This traid is evident from 
Figure 2.6. The reported values of the standard deviation of log(y) ranges from 0.4 to 0.9 

23 Literature Review of Attenuation Relationships 

There are many strong motion attenuation relationships that have been proposed 
throughout the years. Some of the recent relationships, which address the attenuation of 
peak ground acceleration, are summarized below. These empirical equations have been 
arrived at based on different solution procedures, for different types of tectonic regimes, 
site conditions, different distance definitions, and ranges of magnitudes. 

23.1 Boore et oL (1997) 

Boore et al. (1997) studied the peak ground acceleration for shallow earthquakes 
from western North America. The ground motion data consisted of horizontal PGAs of 
271 records obtained from 19 earthquakes. The average of two horizontal components of 
acceleration was used in the regression analysis. The distance range was from 0 to 118 
km and the magnitudes of the recorded events varied from 5.2 to 7.7. Moment magnitude 
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(A/ff) was used as a measure of earthquake size and the closest distance from the station 


to a point on the earth’s surface, which lies directly above the rupture {vjb) was used as 
the distance measure. 


A cutoff distance was introduced in the analysis to account for the triggering of 
instruments due to local site effects at a particular station. The cut off distance was set to 
least of the following: 

a) Distance to the first digitized record triggered by S wave; 

b) Distance to closest non-digitized recording; 

c) Distance to an operational non-triggered instrument. 

Two types of faulting were considered, namely strike-slip and reverse-slip faults. 
The site classification was based on shear wave velocity up to 30m depth. A two-stage 
maximum likelihood method was adopted. The attenuation relationship for the peak 
ground acceleration (g) is given by 


In r = Z), -h 0.527(M - 6) - 0.778 In r - 0.3711n-^ , 
‘ ^ 1396 

where 


( 2 . 6 ) 


- 0.3 1 3 for strike - slip faults 
Zjj = -I - 0. 1 17 for reverse - slip faults 

- 0.242 if mechanism is not specified 

and Fs is the average shear wave velocity which depends on the site class. The standard 
deviation of the predicted acceleration is given as criny= 0.520. The relationship was 
found to give best results when used to predict motions at distances less than 80 km and 
magnitude range 5.5 to 7.5. 


23.2 Campbell (1997) 

Campbell (1997) developed attenuation relationship for mean PGA with majority 
of the earthquakes recorded in California. Some events from Middle East and South 
America were also used in the study. A total of 645 recordings from 47 earthquakes were 
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used for arriving at the attenuation relationship. Moment magnitude {Mw) was used as the 

measure of the earthquake size. The range of magnitudes considered in the study was 4.7- 
8.0. Closest distance from the station to seismogenic rupture {Rsei^ (Figure 2.1) was used 
as the distance measure and its range was 3.0-60.0 km. This distance was used under the 
assumption that fault rupture within softer sediments and within the upper 2 to 4 km of 
the fault zone is primarily non-seismogenic and does not contribute to recorded ground 
motions at periods of engineering interest. 

Two types of faulting, namely strike-slip and reverse-thrust were used to 
characterise the source mechanism. All the earthquakes considered were shallow crustal 
earthquakes (focal depth <25 km). The sites were classified into alluvium or firm soil, 
soft rock, and hard rock and two variables were introduced in the model to represent the 
site conditions. The recommended relationship for PGA (g) is given by 

]n{AH) = -3.512 + 0.904M - 1.328 In 

+ [l . 1 25 - 0. 1 12 \n.{RsEis)- 0.0957M]f 

+ [0.440 - 0.1711n(i?5'£j5)]5'5/{ (2.7) 

+ [0.405 - 0.222 ^(RsEIsi^HR + ^ 

where 

0 .0 for strike - slip faulting 

F = l 0.5 for normal faulting , and 

1.0 for reverse, thrust, reverse - oblique, and thrust - oblique faulting 

Shr = Ser = 0 for alluvium or firm soil; SffR = 0 and -Ss/j = 1 for soft rock; and Se[r= I 

and SsR = 0 for hard rock. 

The regression analysis was done in stages. In the first stage, all selected 
recordings were used to determine the regression coefficients using unweighted, 
generalized nonlinear least-square method. A distance threshold was calculated based on 
magnitude, distance, style of faulting, site conditions, and trigger level of strong motion 
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recorders. In the second stage, the records that did not meet the calculated distance 

. threshold were removed from the set and the regression was repeated. 

Two estimates of standard deviation in the prediction were obtained. The first 
estimate, based on the amplitude of PGA is given by 
0.55, Aff < 0.068g 

cr = < 0.173 -0.140 ln(.%), 0.068g < <0.21g (2.8) 

0.39, Afj>02lg 


The second estimate of cr is based on magnitude Af and is given by the expression 


0.889- 0.069 IM, M <1A 
0.38, M>7.4 


The second estimate for standard deviation was found to be statistically more robust 


2.4.3 Sadigheta/. (1997) 

Sadigh et al. (1997) proposed attenuation relationships for peak ground 
acceleration based on Californian earthquakes. Two major source mechanisms, namely 
strike-slip and reverse faulting were used to classify the records. The relationships were 
developed for rock and deep firm soil deposits. Rock sites are defined as those with 
bedrock less than Im from the surface. Earthquake size was characterized by moment 
magnitude (Af») and distance was defined as the closest distance to rupture surface. The 
distance range of the records used in the analysis was 2.5-261 km and the magnitude of 
the events varied from 3. 8-7 .4. The regression coefficients for mean PGA was calculated 
from 121 recordings from 40 earthquakes. The attenuation model of PGA for rock sites is 
given by 

ln(y) = 6, + b^M + b,{%.5M)^ + b, ln(r„^ + + b, ln(r„^ + 2), (2.10) 

where 

b, = -0.624, b^ = 1 .0, by = 0.0, b, = -2.10, by = 1 .29649, b^ = 0.25 and bj = 0.0 for M < 6.5, 


and 
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b, = -1 .274, ^2 = 1 . 1, 63 = 0.0, b, = -2. 1 0, 65 = 0.4845 1, b^ = 0.524 and b, = 0.0 for M > 6.5, 


with standard deviation given by 

_ ri.39-0.14MforM<7.21 
~ 1 0.38 for M> 7.21 ' 

For deep soil sites the attenuation relationship is given by 

\n{y) = b^+b,M -b^\xi{r^^ + b/^^)+b^+b,{%.5-Mf\ ( 2 . 12 ) 

where 

b\ = -2.17 for strike-slip, -1.92 for reverse and thrust earthquakes; hi = 1 .0; = 1.70 

f2.1863forM<6.5 ('0.32for M < 6.5 

= |o.3825 for M > 6.5 ' "= = |o.5882 for M > 5,5 • ^A, 

with the standard deviation given by 

crinF =1.52-0.16M. (2.13) 


2.3.4 Spudich et aL (1997) 

Extensional regimes are the regions in which the lithosphere is expanding aerially. 
Spudich et al. (1997) developed an attenuation relationship for PGA from worldwide data 
of 129 recordings from 20 earthquakes in the extensional tectonic regimes. Two site 
classifications, namely rock and soil were used. The magnitude of the earthquakes varied 
from 5.1 to 6.9 and the distance (rjb) of recordings ranged from 0 to 102 km. Relationship 
for peak ground acceleration is given by 

logic y = 0.156 + 0.229(M - 6 ) - 0.9451ogio R + 0.077T, (2.14) 

where R =-\lf'jb +5.57^ , and T = 

The standard deviation of the observations is oiogjQ y = 0.216 , which is equivalent to a 
value of 0.497 in terms of natural logarithms. 

2.3.5 Atkinson and Boore (1997b) 

Due to low seismicity rates in stable continental regions, there are very few strong 
motion records available for this tectonic regime. As a result, attenuation relationships are 


0 for rock 

1 for soil 
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obtained using numerically simulated ground motions instead of recorded ground 

motions. 

Based on the numerical simulations for Eastern North America (ENA), Atkinson 
and Boore (1997b) developed an attenuation model for peak ground acceleration. The 
regression is based on thousands of records from hundreds of simulated ENA earthquakes 
in the magnitude range from 3 to 7 recorded at distances (rjb) from 10 to 1000km. The 
digital data used are from hard rock sites. No distinction of focal mechanism type was 
made in the study. The attenuation relationship for PGA (g) is given by 

ln(y) = 1 .841 + 0.686(M - 6) - 0.123(M - 6)^ - Ini? - 0.003 1 li? . (2.15) 

The standard deviation of the observation is found to be equal to 0.62. 

23.6 Youngs et aL (1997) 

Youngs et al. (1997) developed an attenuation relationship for subduction zone 
interface and intraslab earthquakes using the data from Alaska, Chile, Cascadia, Japan, 
Mexico, Peru and Solomon Islands. The rate of attenuation of the ground acceleration is 
less in subduction zone earthquakes as compared to shallow focus earthquakes. The 
attenuation models indicate that for large events at large distances, the peak motions from 
subduction zone earthquakes will be larger than those predicted using attenuation 
relationship for shallow crustal earthquakes. The difference increases as the size of 
earthquake increases. Most subduction zone events are recorded at large distances 
because the events tend to be deep and off shore m most of the cases. In the analysis, the 
events that occurred at a depth less than 50km were considered to be interface 
earthquakes. 

The data used for the regression analysis consisted of geometric mean of two 
horizontal components of 481 recordings from 163 earthquakes. The moment magnitudes 
{Mfv) of the events varied from 5 to 8.4 and the closest distance to rupture surface ranged 
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from 8.5 to 550 km. The regression of the sample was carried out using random effects 

model [Abrahamson and Youngs, 1992]. The random .effects model is a maximum 

likelihood method that accounts for correlations in the data recorded by a single 

earthquake. The attenuation relationship for rock sites is given by 

InO) = 0.2418 + 1 .414M - 2.552 ln(r^^ + 1.781 ) + 0.0060H + 0.3846Zj , (2.16) 

where H is the depth of focus in km and It is the term denoting the source type. It takes a 
value equal to 0 for interface 1 for intraface earthquakes, 
the standard deviation of the predicted motion given by 

ain=1.45-0.LW. (2.17) 

For soil sites the PGA (g) is given by 

ln(y) = -0.6687 + L438 jW- 2.329 ln(r^^ +1.097e®-^^^^) (2.18) 

+ 0.00648/7 + 0.364322’ 

and crin=1.45-0.1M (2.19) 

2.3.7 Crouse (1991) 

Based on the recordings of peak ground acceleration from 697 components 
Crouse (1991) calculated an attenuation relationship for subduction zone earthquakes in 
Cascadia region. The magnitude of the earthquakes varied from 4.8 to 8.2 and the 
distance to center of energy release varied from 8-866 km. The distance to center of 
energy release was assumed to be equal to hypocentral distance for events with moment 
magnitude {Mw) less than 7.5. For larger events, the closest distance to centroid of fault 
plane defined by after shocks was taken as the distance to energy release. The PGA (gals) 
is given by 

In(PG^) = 6.36 + 1.76M - 2.73 ln(i? + 1.58e®-^°^^ ) + 0.00916/7 (2.20) 
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with the standard deviation (crjjj) of the observations equal to 0.773. It was observed 

during the analysis that introduction of a quadratic term of magnitude in the attenuation 

relationship did not have any statistical significance. The solution of the regression 

equation was done using a generalised nonlinear regression code. 

2.3.8 Ambraseys (1995) 

Based on the data collected from earthquakes in Eurasia and Middle East, Ambraseys 
(1995) developed an attenuation relationship for shallow crustal earthquakes. The 
database consisted of 1,667 triaxial records generated during 865 earthquakes. The 
moment magnitude of majority of earthquakes ranged from 2 to 7.3 and the recorded 
distances ranged from 0 to 260 km. For larger events in the database {Ms > 6.0) and for 
earthquakes for which results from special focal studies are available, d=rjb was used as 
the distance measure. For smaller events epicentral distance {d=lQ was used as the 
distance measure. When the events with magnitude > 4.0 was considered for the analysis 
(828 records from 332 earthquakes) the following empirical equation for peak ground 
acceleration (g) was obtained. 

log(y) = -1.06 + 0.245^5 -0.00045r-1.0161og(r), (2.21) 

where r = and h is the depth of focus of the earthquake. 

The standard deviation of the prediction is given by criogjo ” 0.254, which is equivalent 

to CTin =0.585. 

Two-stage regression procedure was adopted for the analysis of the data. It was 
also observed that the predicted attenuation curves using one stage regression and two- 
stage regression do not differ very much. One-stage regression values were slightly 
higher at distances less than 10 km and slightly lower at distances greater than 100 km. 
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2.3.9 Tento et aL (1992) 

The data from Italian earthquakes until 1991 was used by Tento et al. (1992) to 

develop the attenuation relationship for PGA. A two-stage regression anah^is was carried 

out. The local magnitude was used in the study with hypocentral distance as the distance 

measure. The resulting equation is as follows 

ln(y) = 4.73 + 0-52M - Ini? - 0.002 16i? , (2.22) 

where R = +H^ and H is the depth of focus. 

The standard deviation ( (Tin ) the model is equal to 0.67. 

2 J.IO Fukushima and Tanaka (1990) 

Fukushima and Tanaka (1990) developed an attenuation relationship for near 
source region in Japan based on 1,372 horizontal components of PGA from 28 
earthquakes in Japan and 15 earthquakes from other countries. The estimated ground 
motion equation is given as 

logio y = 1 .3 + QAlMs - logio(i? + 0.032 x lO®-'^^^ ) - 0.0034i? , (2.23) 

where R is the shortest distance between site and fault rupture {Rn^. 

The standard deviation of the model was o'logjQ = 0.391 i.e., ctjjj = 0.90 . 

23.11 Caillot and Bard (1993) 

Data collected from Italian sites were used to establish empirical relations 
between spectral acceleration, magnitude, hypocentral distance and site conditions. Total 
of 115 records from earthquakes of magnitude between 3.2 and 6.8 with epicentral 
distance less than 60 km and focal depth less than 30 km were used. The records 
belonged to three different site categories namely, rock (So), thin alluvium with thickness 


less than 20m (Si), and deep alluvium (S 2 ). The form of the equation is given by 
ln(y) = -3.692 + 0.507M - OASSH^hypo) + 0-2715' , 


(2.24) 
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where 5 = 0 for site category So and 1 for site category Si. Since the data coming under 

site category S 2 was too heterogeneous, it was not considered in the study. 

For the data considered in the study, regression analysis was done using two 

multi-linear regression methods; (i) two stage regression analysis which de-couples the 

spatial dependence of acceleration from that of magnitude, and (ii) weighted multi-linear 

regression procedure where weighing is aimed at correcting the non-uniform magnitude- 

distance distribution for each site-specific data set. It was observed that the two-stage 

regression procedure does not reduce the initial variance (stage 1) of the data 

< 0). It was due to the fact that the two successive regressions correspond to a 
local non-absolute minimum and not the numima of the whole surface. 

13.11 Anders eta/. (1990) 

Andres et al (1990) developed an attenuation relationship based on the 87 strong 
motion recordings from 56 intraplate earthquakes that occurred in North America, 
Europe, China and Australia. A two-stage regression was used. Hypocentral distance was 
taken as the distance parameter and the surface wave magnitude (AG) was used as the 
measure of size. The largest horizontal component of the horizontal PGA was used for 
the analysis. The relationship is given by 

In ^ = -1.471 + 0.849M - 0.00418i? + lnG(i?,i?o) (2-25) 

R~^ forR<R^ 

with o-jn =0.984. Here, G(i?,i?o) = ] -The recommended 

Rq 1^1 for R ^ Rq 

value for Ro is 100km. 


• • • 
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Table 2.1: Comparison between the parameters obtained by one-stage maximum 
likelihood method, the two-stage method and the ordinary least squares 
method applied to Fukushima and Tanaka ’s data for the model 
log(>^) = a + bM + c logi? [Joyner and Boore, 1993] 


Parameter 

One-Stage 

Maximum-likelihood 

Two-stage 

Maxim um-likelihood 

Ordinary 

Least-squares 

a 

2.20 

2.20 

2.22 

b 

0.44 

0.44 

0.26 

c 

-1.73 

-1.76 

-1.19 


0.26 

0.26 

- 



0.21 

- 

& 

0.33 

0.34 

0.33 
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Vertical faults 



(a) 


Dipping faults 


rjb=0 



Figure 2.1: Different definitions of distances used in the attenuation relationships 
[Abrahamson and Shedlock, 1997]. 
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Ml: Hypocentral Distance; 

M2: Epicentral Distance; 

M3: Distance to Zone of energy release; 

M4: Closest Distance to rupture; 

M5: Closest Distance to surface projection of Rupture 

Figure 2.2: Different definitions of distances used in the attenuation relationships 
[Campbell, 1985]. 
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MOMENT MAGNITUDE 


Figure 2.3; Relationship Between moment magnitude and various magnitude scales 
[Kramer, 1996]. 
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DISTANCE COEFFICIENT 


Figure 2.4: Histogram showing distribution of distance coefficients for individual 
earthquakes and average value along with coefficients obtained by non- 
stratified (one step) and two-step analyses [Fukushima and Tanaka, 1990]. 
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(a) 

EAfTHQUAKE PEAK GROUND ACCELERATION 



Comparison of predictive curves for peak acceleration from two-stage and 
one-stage regression, (a) Predicted equations without considering focal 
depth, (b) Predicted equations when focal depth is considered in the analysis 
[Ambraseys, 1995]. 



standard Deviation (a) 
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4.5 5 5.5 6 6.5 7 7.5 8 8.5 

Magnitude (M) 


Figure 2.6: Variation of Standard Deviation ( ) of Peak Ground Acceleration with 

respect to magnitude for different attenuation relationships. 


Chapter 3 

Strong Motion Data from Indian Earthquakes 


3.0 Introduction 

The Himalayan region in northem-India is a highly seismically achve region. 
Because of the economic as well as operational constraints, only a small area in northem- 
India could be covered with Strong Motion Accelerographs (SMA’s). Stmctural 
Response Recorders (SRR) are relatively inexpensive instruments, which are used to 
record the spectral accelerations during an earthquake. Due to low cost, SRRs have been 
installed covering a wider region. This chapter compares the characteristics shown by the 
recordings from SMA’s and SRR’s. A methodology is also developed to calculate Peak 
Ground Acceleration (PGA) from Resultant Spectral Accelerations (RSA’s) recorded by 
SRR’s. 

3.1 Seismo-Tectonics of Himalayan Region 

The northeastern region of India is an active seismic zone. An earthquake of 
magnitude 8.7 has occurred in the Shillong region in 1897 and in the Mishmi region of 
Arunachal Pradesh in 1950. Shillong plateau is bounded by the main boundary thrust in 
the north-west, and the Mishmi and the Lohit thrusts towards the north-east (Figure 3.1). 
The other two prominent tectonic features forming the boundary of the Shillong plateau 
are the Dhubri fault towards the west and the Dauki fault towards the south. The Dauki 
fault merges with the north-east trending Haflong-Disang thmst towards the east. A 
thrust type mechanism was suggested for the Dauki fault in a recent study [Qiandra and 
Das, 1992]. There is also report of a subduction process along Burmese arc, which lies 
near the Shillong region [Mukhopadhyay, 1992]. 

The Bihar region consists of a tectonic basin in the east-west direction along the 
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southern margin of the Himalayan fold belt [Dasgupta, 1993]. The formation of this 
basin is attributed fo the upheaval of the Himalayas along with the down warp of Indian 
lithosphere. The central part of this basin is called the Gangetic foredeep. The 1988 
Bihar-Nepal A/6.8 earthquake originated below the Gangetic foredeep [Dasgupta, 1993]. 

In the Kangra region, the major tectonic features are the Main Boundary Thrust 
(MBT) and the Main Central Thrust (MCT) (Figure 3.1). 

3.2 Instrumentation and Availability of Data 

The Himalayan region of India has been instrumented in the recent times using 
the Strong Motion Accelerographs (SMA) and the Structural Response Recorders (SRR). 
Strong Motion Accelerographs are three-component instruments installed near a 
potential seismo-tectonic source; these provide time history of ground acceleration in 
two horizontal and one vertical directions. In all, there are five major arrays of SMAs in 
the Himalayan region [Chandrasekaran et al, 1996]. Figure 3.2 shows the relative 
positions of the strong motion arrays installed in different parts of the Himalaw region. 
Under the Indian National Strong Motion Network (INSMN) Program, Department of 
Earthquake Engineering (DEQ), Roorkee has installed more than three hundred SRRs in 
different parts of India. The multiple Structural Response Recorder (SRR) consists of six 
recording instmments, in the form of conical pendulums, which measure the horizontal 
component of motion (Figure 3.3a). These seismoscopes are relatively inexpensive 
instruments. Their natural periods are adjusted to 0.40sec, 0.75sec and 1.25sec and the 
damping provided is 5% and 10% of critical damping. The output from a SRR provides 
ordinates of Resultant Acceleration Spectrum (RAS) at natural periods 0.4sec, 0.75sec 
and 1.25 sec for 5% and 10% damping. A typical SRR recording is shown in Figure 


3.3(b). 
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3.2.1 Shillong Array in Meghalaya 

This array consists of 45 strong motion accelerographs located between 24°N to 
26°N and 91°E to 94“E (Figure 3.4a). About 110 SRR stations are also situated between 
23°N to 28°N and 90°E to 95°E. Strong motion accelerograms recorded by the Shillong 
strong motion array for five earthquakes are available, viz., 10 September 1986 
Meghalaya (mi=5.2) earthquake (12 records), 18 May 1987 Northeast {Ms=5.9) 
earthquake (14 records), 6 Feb 1988 northeast (/ni=5.8) earthquake (18 records), 6 
August 1988 (Ms=7.3) earthquake along the Indo-Burma border (33 records), and the 9 
January 1990 Northeast (m 6 = 6 . 1 ) earthquake (14 records). The 6 August 1988 
earthquake was also recorded by about 120 SRRs installed in the region. Figures 3 . 5 (a) 
to 3.5(e) show the locations of strong motion recorders and SRRs, which have recorded 
these earthquakes along with the epicenter locations as per different agencies. 

3.2.2 Bihar Array in Bihar 

This array (Figure 3.4b) consists of five digital accelerographs (Chandra et al, 
1996) and about 47 SRRs. The 21 August 1988 Bihar-Nepal earthquake was recorded by 
41 SRRs installed in the region (Figure 3.5f). 

3.2.3 Uttar Pradesh Array in North-Western Uttar Pradesh 

This array contains 40 SMA stations between 29®N ’to 31°N and 77.5°E to 81°E 
and about 70 SRR stations between 28“N to 32°N and 76.5°E to 80.5°E (Figure 3.4c). 
The 20 October 1991 Uttarkashi (M=7.0) and 29 March 1999 Chamoli (M 5 = 6 . 6 ) 
earthquakes occurred in the seismic gap lying between the rupture zones of the Great 
Kangra (1905) and the Bihar-Nepal (1934) earthquakes. Uttarkashi earthquake was 
recorded by 13 SMA stations of the Uttar Pradesh array (Figure 3.5g). 1 1 1 SRR stations 
including those situated at far off places like New Delhi also produced records of this 
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event. The Chamoli earthquake was recorded by 11 SMA stations and 16 SRR stations 
(Figure 3.5h). 

3.2.4 Kangra Array in Himachal Pradesh 

This array (Figure 3.4d) contains 50 SMAs located between 31'^ to 33'Tl and 
75°E to 78°E and about 65 SRR stations between 30.5°N to 33.5^^ and 75.5'’E to 78.5°E. 
The 26 April 1986 earthquake {mb=5.5) was recorded by 9 SMAs [Chandrasekaran and 
Das, 1992] and 17 SRRs (Figure 3.5i). 

3.2.5 Arunachal Pradesh Array in Arunachal Pradesh 

The Arunachal Pradesh array consists of 20 digital accelerographs (Chandra et. 
al, 1996). 

The earthquakes of which the data is analyzed in this study are listed in Table 
3.1(a). The numbers of strong motion recordings obtained along with the range of 
epicentral distances are listed in Table 3.1(b). The peak parameters like maximum 
acceleration, velocity, and displacement for these earthquakes are given in Table Al and 
SMA station coordinates are given in Table A2 of Appendix A. The values of spectral 
acceleration recordings by SRRs are listed in Table B1 and the SRR station coordinates 
are given in Table B2 of Appendix B. 

3.3 Earthquake Parameters 

Wherever available, locations of epicenters proposed by United States Geological 
Survey (USGS) and Indian Meteorological Department (IMD) for different events are 
shown in Figures 3.5(a) to 3.5(i). For some of these events, epicenter location obtained 
by Chandrasekaran and Das (1990) using the strong motion records are also shown; these 
are denoted as DEQ. Epicenters by different agencies differ significantly and sometimes 
the difference between these can be up to lOOkms (e.g., Figure 3.5c and 3.5d). Table 
3.1(a) shows the earthquake magnitudes calculated by different agencies. As the 
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magnitude calculations for all the earthquakes considered in this study are available from 
USGS, the magnitude as reported by USGS has been followed. 

For all the earthquakes, the correlation coefficient betv.een the average PGA in 
the two horizontal directions and the epicentral distance was calculated (Table 3.2). The 
epicenter that gave the best correlation coefficient has been considered in the present 
study. 

For most earthquakes, the causative fault details are not available in the literature. 
Similarly, there is a significant uncertainty in the focal depth. Hence, the epicentra. 
distance is used as the control distance to the fault. 

3.3.1 Calculation of Epicentral Distance 

From known coordinates of epicenter and recording station, the epicentrai 
distance should normally be calculated by assuming shape of the earth as spheroid. 
However, assuming shape of the earth as a sphere introduces an error less than 0.21% for 
a distance of 400 km. The error will be less for smaller distances. Hence, the spherical 
triangle formula was used for the calculation of epicentral distances. 

The spherical angle ZMN between two points M and N on the surface of earth, 
with latitude and longitude {La^,Lojfi) and (La„,LOn), respectively, is given by (e.g., 
Abramowitz and Stegun, 1972): 

cos(ZM/V) = cos(90 - La jn) cos(90 - ) + sin(90 - La „) sin(90 - Laj„ ) cos(Lo„ -Lo„) 

(3.1) 

The mean radius at a point on surface of the earth for the mean latitude <(>^ is given b>' 
[e.g., Clarke, 1950] 


a\ 




J/2 


(3-2) 
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where the values a— 6,377,304 and A— 6,356,103 are used in India by Geological Survey 
of India (Everest’s values) and =l-(h/u)^. The mean latitude of stations ({zij in 
Shillong array was found to be 25°9' and this value was used in the calculation of as 
6346.38km. It was observed that the change in Rm due to the change in mean latitude 
from 25°9' to 32° (Kangra array) was rather low (<0.1%) and hence the same value of 
was used for all the events. For the sphere of radius 6,346.38km, the epicentral distance 
in km is given by: 

R = 63463SxZMN^ (3.3) 

3.4 Analysis of SRR Data 

As seen in Table 3.1(b), for several earthquakes considerably more number of 
SRR recordings are available than the SMA records. Therefore, any study on ground 
motion attenuation must consider both SMA and SRR data. A systematic methodology 
has been used in this study to estimate PGA values from SRRs, and thus the database for 
PGA values has been considerably enhanced. However, SRR is a relatively simple 
instrument and it is necessary to first verify the quality of SRR data. Rest of this chapter 
is concerned with the following objectives: (a) To assess the degree of uniformity 
between the SRR data, (b) To examine the similarity between information from SMA 
and SRR data, and (c) To develop a methodology for the computation of PGA from SRR 
data. 

3.4.1 Resultant Response Spectrum 

The SRR produces records in a two-dimensional reference firame, i.e., any point 
on the SRR record results from the two-dimensional input in the horizontal plane. 
However, the SMA records are unidirectional ones; the input motion is recorded along 
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mutually perpendicular directions. Thus, before the records from SMA and SRR are 
compared, the SMA records have to be reduced to represent the two-dimensional input. 

The time history record of ground acceleration along longitudinal direction was 
used as an input to calculate the response acceleration of a Single Degree Of Freedom 
(SDOF) oscillator with natural period T and damping of critical. This was done with 
the help of program SPECTR [NICEE, 1998]. This procedure was repeated using the 
ground acceleration along transverse direction. The resultant of these response 
accelerations acting along mutually perpendicular directions was calculated by taking the 
vector sum of response accelerations at each time instant. This procedure is repeated 
covering the complete time history. The maximum value of resultant response 
acceleration obtained during the entire time history is taken as the Resultant Spectral 
Acceleration (RSA) for natural period T and damping (^% of critical. Figure 3.6 shows 
the schematic representation for the calculation of the RSA from SMA data that can be 
directly compared with SRR output. 

3.4.2 SRR Data versus SMA Data 

Since SRR record provides maximum response acceleration in any direction, it 
gives ordinates of the resultant response spectrum. These ordinates are provided at 
periods 0.4sec, 0.75sec and 1.25sec for 5% and 10% damping, i.e., 0.4sec and 5% 

damping (represented as SAq 4 ^§o/^), at period 0.75sec and 5% damping (‘?-4o.75^5o/„) etc. 

As the SMA and SRR are seldom located at the same location, the direct 
comparison of SRR and SMA data is not possible. However, if the attenuation 
characteristics for the region are similar for all the recording stations, the statistical ratios 
of recordings should show similar trends. 

Nine different ratios of SRR recordings are studied and compared with those 
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damping ratio, for instance between periods 0.4sec and 0.75sec for 5% damping (denoted 

as ( 54^4 / 'The remaining three ratios studied are between the RSAs 

for the same period but for different damping values, e.g., between period 0.4sec for 5% 
damping and period 0.4sec for 10% damping (denoted as (5^5o/o/‘S'4Vo)r=0 4 )’ 

3.43 Removal of Outliers 

The various ratios for SRRs and SMAs as discussed above were plotted to form 
histograms. These histograms can be divided into two categories: (1) Ratios of records 
with same damping ratio but with different natural periods, and (2) Ratios of records 
with same natural period but with different damping ratios. The histogram for a 
particular ratio for SMA and SRR are given side by side for ease of comparison. (Figure 
3.7a to 3.71). It can be observed from these histograms that: 

a. The shape of the frequency distribution followed by SRRs as well as SMAs is 
similar. 

b. For both SRRs and SMAs, the peaks of the distribution lie almost near the same 
value. 

c. The shape of the distribution functions are either symmetrical or skewed to the 
left. It can be assumed that these ratios follow either normal or lognormal 
distribution. 

d. The distribution follows a well-behaved form for the ratios of RSA corresponding 
to (s 4 _ 4 /s 4 js\ and [sAojsf dispersion in the histogram 

increases for the values of This shows that, the ratios of RSAs 

with closely spaced natural periods show less scatter compared to the ones with 
much larger difference in natural frequencies. 
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e. For some SRR stations the ratio of response with 5% damping to that with 10% 

. damping is less than 1.0. However, the response of a 5% damped system should 
be greater than that with 10% damping. Similarly, for some SRR stations the ratio 
appears too high. 

The unusually high or low values might have resulted from some erroneous 
records and poor calibration of the SRR instruments. Hence, it is desirable to remove 
these records or outliers from the data sets before they are used in further analysis. 

Figure 3.8 gives the frequency distribution of the average value of 

(‘5'-^0.4sec,$'=5% /'^■^0.4sec,$-=10% )’ (‘^■^0.75 sec, ^=5% / ‘^0.75sec,5'=10% )> 

(■^■^1.25sec,f=5%/‘5'-^1.25sec,$-=10% ) SMAs and SRRs. For August 6, October 20 and 

August 21 events the mean value of the above frequency distribution for SRRs was 
found to be 1.5,1.49 and 1.4 with standard deviations equal to 0.4, 0.37 and 0.3. In the 
case of SMAs, for August 6, October 20 events, these mean values were 1.3 and 1,24 
with standard deviations 0.06 and 0.07 respectively. It can be seen that, the mean values 
obtained from SMAs are less compared to that of SRR and SD of the ratios from SRRs is 
very high compared to that of SMAs. One of the reasons for this xmusual difference in 
the values may be due to inclusion of erroneous observations. In order to remove these 
points, the following strategy was adapted. 

The average of mean + a values of these events was found to be equal to 1.83. 
So the upper cut-off point was fixed at 1.85 (Ucut) and the lower limit at 1.0 (Z.cu/)- Based 
on this criteria, about 18 (l6Ucut+2Laa) observations were removed from August 6, 15 
(15l7cu,) from October 20 and 3 (lUaa+lLcu) from August 21 event. Removal of points 
for which the ratio do not lie between 1.0 & 1.85 had the following effects; 

• The change in mean value of [sA^yJ pronounced compared to the 
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Ucuh it meant curtailment of the tail of the distribution and this had more pronounced 
effect on mean than for median. 

• The change in various statistical parameters of response ratios between different 
damping values (e.g., SA^jffy^) was more significant compared to that for 

response ratios across different periods (e.g., [sA^JsA^^ ) ) for the same amount of 

damping. This may be due to the fact that the removal criterion is based on the 
response ratio for different damping values. 

The rest of the analysis uses the data remained after the above removal of outliers. 

3.4.4 Observations on Data Analysis 

To gain more confidence on the SRR data, correlation coefficients (CC) between 
the variables considered in the analysis were calculated for SRR as well as SMA data. 
These along with the other statistical coefficients are given in Table 3.3(a) to Table 
3.3(e). For the events considered in this study, the following general observations can be 
made; 

1. The mean and median values of RSA ratios from SRR data are fairly close to those 
from SMA data (e.g., Tables 3.3a to 3.3e). Generally, the median values from two 
data sets compare better than the mean values. The average difference between 
median values of all the ratios is about 16% and the difference ranges from 1% to 
68% whereas the average difference between mean values is less than 20%; and it 
ranges from 3% to 45%. 

2. The mean ratios of RSAs from SRRs for different damping corresponding to a 

constant period )r generally higher compared to the corresponding 

ratios obtained from SMA records (average difference of 19%). 
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3. The median values of the ratios of RSAs from SRRs based on different damping 

values good agreement with that of SMAs (average difference 

of 5%). At the same time, the difference in the median values of ratios of RSAs 
between different natural periods with same damping ) shows 

considerable difference (average difference of 25%). 

4. The standard deviation of RSA ratios for different natural periods but with same 
damping from SRR data compares well with that from the SMA data. However, 
standard deviation for RSA ratios for different damping values but same natural 
period is much higher for the SRR data than that from the SMA data; in fact, the 
value from SRR data may be four to five times that from the SMA data. This is 
perhaps an intrinsic nature of the data. In case of SRR different oscillators record the 
motion for 5% and 10% damping values and there tends to be a large variation in the 
ratio of the two quantities. For SMA, same time history is used to computationally 
obtain the response for 5% and 10% damping ratios and thus there is only a limited 
variation in the two quantities. 

5. The correlation factors of the ratios of RSAs from SMAs and SRRs show good 
agreement with each other for all the earthquakes. However, for 29 March Chamoli 

earthquake, the correlation coefficients corresponding to 


(^•^.75/ ^^\. 25 \= 5 % 


for SRRs show large deviation from the corresponding values 


for SMA. This may be due to poor calibration of the instrument with period 1.25sec 
and 5% damping. 

From the above observations, it can be concluded that the characteristics shown 
by SMAs and SRRs are similar in nature. The similarity between the recordings of 
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SMAs and SRRs has also been studied with the help of statistical tools in the following 
section 

3.4.5 Statistical Analysis of Data 

Many of the statistical tests available need certain assumptions to be made on the 
nature of the distribution followed by the variables. From the plots of response ratio 
histograms, it can be observed that their shape resembles a normal distribution. The 
frequency distribution of the ratios for RSAs of SRRs obtained in the August 6, 1988 
earthquake was tested to find out validity of this observation. A test which can be used to 

check the validity of the model is ti?" test for goodness of fit [Kreyszig, 1979]. The 


{^^test was conducted for the sample under consideration and the assumption of 

distribution being nonnal was found accqptable with 5% level of significance. 

When the underlying populations follow normal distribution, the following test 
statistic can be used to test the fact that the response ratio values obtained from SRRs 
and SMAs have originated from the same population [Miller, F., 1977]. 


1 = 


(^1 




(3.4) 


where, 3c; , < 7 ; are the mean and standard deviation of the samples obtained from SRRs 
and X 2 ,ct 2 those of SMAs and 2^^® corresponding number of observations. ^J’is 

the difference in the population means ///, pti which will be equal to zero if the 
population means are same. The assumption that the samples of SRRs and SMAs do not 
come from the same parent population can be considered acceptable if test statistic (z) do 
not lie in the range +1.96 to -1.96 (with level of significance, a = 5%). The results of 

statistical analysis are given in table 3.3 (a) to 3.3 (e). The ratios j 
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the hypothesis about the parent population of the two samples being same. This 
reinforces the assumption that the SRRs as well as SMAs are generated from the same 

parent population. The hypothesis is generally not satisfied for for the 

data considered, except for those from the 26 April 1986 Kangra earthquake. This may 
be because for all the other events the standard deviation for the response ratios of SMAs 

with same period and different damping are very low (up to 1/5^*) as compared to that 
from the SRR data and the possible reasons for the same have been discussed earlier. 

Based on the results of this comparison, one gains more confidence on the SRR 
data. However, one must avoid using SRR data for study of response ratio for different 
damping values. 

3.4.6 Spectral Acceleration versus Distance 

The variation of RSA with distance is studied by calculating the correlation 
coefficient (CC) betv.een the observed RSA (from the SMA and well as the SRR data) 
and the epicentral distance (Table 3.4). This parameter can also be used to bring out the 
similarities/differences between SMA and SRR database. The following observations 
can be made: 

1. The CCs calculated for SRRs and SMAs are quite comparable. The average value of 
CC is around -0.5, 

2. For August 6 , 1988 Shillong earthquake, the average CC for SRRs (-0.523) is found 
to be higher compared to that for SMA (-0.484). 

3. For Uttarkashi and Chamoli earthquakes, the CCs for RSAs of SRRs (-0.397 & - 
0.478) are showing lower correlation when compared with the corresponding values 
for SMAs (-0.613 & -0.549, average difference being 30% and ranges from 15 to 

SJo 


55 %). 
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4. In the case of Kangra earthquake, there is high correlation between the SMA spectral 


values and distance (CC«-0.9). But, SRR records do not show this trend (CC«-0.4). 
The difference between the observed correlation factors for SRRs and SMAs is 
maximum in this case (average difference of 55%). It is also observed that, for this 
event the CC corresponding to the response from the SRRs with r=0.4 and ^=5% is 
very low (a-0.1) but the CC calculated from rest of the SRRs give almost uniform 
value (»-0.45).. 

5. In northern Himalaya, the RSAs from SRRs are less correlated with epicentral 
distance compared to that from SMAs. 

For most of the earthquakes CCs of SRRs are close to CCs of SMAs with respect to 
epicentral distance and this nature also strengthens the observation that SRRs and SMAs 
show similar properties. 

3.5 PGA from SRR Data 

Spectrum Amplification Factors (SAF) are commonly used to develop response 
spectrum from known PGA values. Newmark and Hall (1982) have proposed values of 
amplification factors for various amounts of damping (Table 3.6) and methodology of 
calculation of response spectrum bounds employing these values. This was based on the 
fact that response spectrum over certain frequency ranges is related to the peak value of 
ground acceleration, velocity and displacement. For calculating PGA from resultant 
spectral values a s imilar concept can be used. But, Chandak and Jain (1994) have 
reported that the amplification factors are significantly different from that reported in the 
study by Newmark and Hall (1982). 

A methodology can be developed based on the characteristics of the recordings 
from SMAs and the validity of the procedure can be tested based on the available values 
of actual PGAs for SMAs. The procedure which gives the best prediction of PGAs can 
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be employed in conversion of SRR responses to PGA values. Before amplification factor 
is used as the conversion parameter, the-variation of amplification factors with respect to 
distance and with change in geological characteristics of different regions needs to be 
studied. 

3.5.1 Spectrum Amplification Factor 

The high frequency waves of the ground motions attenuate faster than low 
frequency waves. This may lead to SAF being influenced by distance. In general, 
SAIPGA values are generally expected rise with respect to distance in velocity and 
displacement controlled region. The variations of Amplification Factor (AF) with 
distance along with the logarithmic best-fit lines are shown in Figure 3.10 (a) to 3.10 (d). 
Table 3.5 shows various statistical parameters calculated between RSA and average PGA 
from SMAs. Table 3.7 shows the correlation coefficients of various amplification 
factors with respect to epicentral distance {SAjPGA , D)j- ^ . 

It can be observed from Table 3.7 that CC changes from +0.53 to -0.501 for 
different events and there is less negative correlation between amplification factor and 
epicentral distance. It is seen that the amplification factors corresponding to the same 
earthquake shows an increasing trend with respect to distance for some natural periods 
and decreasing trend for some other values (e.g., Uttarkashi, Chamoli). For some events, 
the trend is consistent for all periods (e.g., Shillong). The relationships for spectral 
acceleration available in the literature predict the relation of amplification factor with 
distance as one with a very weak negative dependence or a strong positive dependence 
(Figure 3.11). The trend predicted by all the relationships is consistent for periods 0.4sec, 
0.75sec and 1.25 sec [Youngs et al, 1997; Atkinson and Boore 1997; Sadigh et al, 
1997; Boore et al, 1997], which contradicts with the observations made for Indian 
earthquakes. The logarithmic best fit for the data of amplification factors versm 
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epicentral distance tends to become flat (weak dependence) as the period and damping 
increases. Thus, from the analysis of correlation factors between {SAlPGyi)j .and D and 

the trend of amplification factors {SAj PGA)j ^ versus epicentral distance Z), no general 

conclusion can be drawn about the distance dependence of amplification factors. The 
records from SRRs cover large distances compared to SMAs. Thus, any distance 
dependant extrapolation of data based only on the observations from SMAs might lead to 
erroneous results. Hence {SAj PGA)j ^ is assumed independent of epicentral distance D 

for the regions considered in this study. 

Since the amplification factor is considered independent of distance, the regional 
variation of the amplification factor can be studied with tie help of parameters given in 
Table 3.5. The calculated parameters include mean ratios of RSA and average PGA at 
different time periods and at different damping ratios, represented as {SAfPGA)j-^, and 

their medians, standard deviations and correlation factors. The following observations 
can be made on these parameters corresponding to different events; 

1 . The mean SAF for these events ranges from 2.45 to 2. 14, 1 .0 to 1 .44 and 0.45 to 0.77 
corresponding to 5% damping for natural periods of 0.4, 0.75 and \J5 sec. The 
corresponding values for 10% damping were found to range between 1.68 to 1.87, 
0.84 to 1.07 and 0.38 to 0.61. 

2. Based on the local geology, the mean values of amplification factors show' significant 
variation from event to event (e.g., AFi.25,io%=0.377 for Kangra and 0.605 for 
Shillong earthquake). Similar trend is also shown for median values. This implies 
that the shape of the normalised response spectrum will be different for each event 
and use of a single amplification factor corresponding to all events for a particular 
^|;riod will be meaningless. So, in the subsequent analysis the amplification factors 
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used in the computation of SRRs are the ones obtained from the SMAs of 
corresponding events. 

3. Generally, median values are lower than the mean values. 

4. The standard deviation of SAP decreases as the natural period and damping ratio 
increases. But, the rate of decrease is different for each earthquake. Since the mean of 
amplification factor also decreases with increasing period, it is observed that the 
coefficient of variation increases with period and damping. 

5. Higher value of damping produces better correlation of amplification factor with 
distance. 

6. For 26 April 1986 Kangra earthquake and March 20, 1999 Chamoli earthquake PGA 
values were found to have very good correlation with RSA (ranges firom 0.82 to 
0.94). This implies that the ratio of RSA/PGA is not being influenced by any external 
elements like geology. The least correlation with PGA was observed for 6 August, 
Shillong earthquake. For this event, good correlation between RSA and PGA (0.88) 
was observed for RSA corresponding to natural frequency of 0.4. But, these were 
found to decrease as the natural period of the system increased (0.6 for r=0.75 and 
0.45 for r=1.25sec). 

Based on the above observations it can be seen that the mean or median values of SAF 
can be used for arriving at PGAs. The use of CCs as weighing functions can also be 
explored. This aspect is discussed in detail in the next section. 

3.5.2 PGAs from Recorded RSAs 

Different alternatives for calculation of PGA from RSA are explored in this 
section. All methods are based on employing SAFs to back-calculate PGAs. Since the 
actual PGA values are available for SMAs, the accuracy of the different procedures can 
be compared. The parameter known as error of estimate was used to compare the 




When this procedure was employed, the error in prediction for different events ranged 
from 9.4 to 17 % (Table 3.8). 

3. 5.2,2 Method Based on Median of Analysis and Weight Functions 

As the correlation coefficient relates the degree of dependence of RSA on PGA. 
this method assigns the RSA values with higher correlation with PGA more weight in the 
calculation of the PGA, i.e., the correlation coefficients are used as weighing function for 
each RSA. The peak ground acceleration is obtained by dividing the sum of obser\'ed 
accelerations by the sum of weights used and written as 

I median 

PGA predicted ~~ ^ • (3-7) 

i=l 
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It is seen that the error ranges from 9.4 to 16.9 % in this procedure (Table 3.8). ' 
reduction in error is considerable in the case of August because compared to others, 
observation with 0.4 sec period had shown more correlation with PGA. But for all ot 
events this method did not result in any significant reduction in error. This is due to 
fact that the CCs for these events are almost equal to unity and hence this method v 
not produce results significantly different from the previous method. 

3.5.2.3 Method Based on Mean of Analysis 

Tthe mean of the {SAfPGAY^ ratios of the available records is used as t 

amplification factor and weighing functions are assumed to be unity. This can be writt 
as 



where, i represent to each response from SRRs. 

The use of this method resulted in significant reduction in error. For some events, tl; 
reduction was upto 4%. The error of prediction raged from 9.0 to 14.0. 

3.5.2.4 Method Based on Mean of Analysis and Weight Functions 

In this method, the mean is used as the amplification factor and the correlatio: 
coefficients between RSA and PGA are used as the weighting functions. The relation i, 
given as 


PGA 


predicted 



(3.9) 


Out of all the methods discussed above, the least values of errors were obtained when 
this procedure was used. The error in prediction varied from 8.9 to 14. 
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A method based on using arbitrary weighing factors 0.90, 0.60 and 0.45 
corresponding to periods 0.4, 0.75 and 1.25sec was also employed. The use of these 
weighing factors resulted in lesser value of error (Table 3.8). The error ranged nom 8.5 
to 13. For August 6 , 1998 event, the values the weighing factors were approximately 
equal to the correlation coefficients, but these weighing functions were physically 
meaningless for all other events. 

Since a best-fit equation will be providing the least error of prediction, an attempt 
was made to compare results from best fit method with the procedures discussed above. 
3.5.2,5 Best-Fit Method 

In this method a best-fit expression is obtained for the sample considering PGA 
as the dependent variable and the various SA as the independent variables. There are 6 

values for the resultant SA from the SRR records, i.e., SAL for different values of 

period and damping. In this study, two such models are considered. 

Model I : 

PGA = +^> 2 ‘ 5 '-^ 0 . 75 , 5 % ■'■^ 3 ‘ 5 '-^[. 25 , 5 % 

■'■^ 4 ‘^-^ 0 . 4 , 10 % "^^ 5 ‘^•^ 0 . 75 , 10 % "^^ 6 ‘^•^ 1 . 25 , 10 % 

Model II : 

PGA = bi 

+^2 ‘^•^0.4,5% ■•■^3 ‘^■^0.75, 5% ■‘■^4‘5''^1.25,5% 

+^ 6 ‘ 5 -^ 0 . 75 , 10 % * 7 ‘^•^[. 25 , 10 % 

For the model I, the values of the coefficients calculated for August 6 1988 event are 
given below; 

bi = 0.0818 , b2 = 2.988, 63 = -6.303 , b^ = 0.33, *5 = -4.365 , be = 8.7965 . 

Theoretically, these values should be similar to AFj6 and should be greater than zero. It 


(3.10) 


(3.11) 
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can be observed from these coefficients that no general trend can be predicted about the 
relation of this numbers with the Resultant Spectral Accelerations. The usage jof these 
coefficients as weighing functions was also discarded because of arbitrarily high values 
corresponding to some accelerations and change in signs which would imply a negative 
effect of a particular RSA on PGA. 

For the second model described in least-square analysis, the coefficients obtained for 
August 6 1988 event are 

bi = 0.0533 , bj = 0.4748, b-^ = 0.5061, b^ = -0.2177, 65 = 0.9587 , 65 = -0.4820 
b-^ =-0.2835. 

In this case, the calculated error was found to be 5.85%. For some events, model -1 
predicted the least error and model -2 for some other events. 

The performance of the solutions of least-square estimates when extrapolated is 
known to be erroneous especially when there is no physical basis for selection of 
parameters (Campbell 1985). Since the range of recordings by SRR are much wider 
compared to that of SMAs and model-2 do not have any physical basis, the use of this 
method for calculation of PGA from SRR values may produce misleading answers. 

It was also observed that the error predicted by this procedure was not very 
different from the one predicted using mean as amplification factors and correlation 
coefficients as weighing functions (difference of 0.6 to 5). So for the back calculation of 
PGA from SRRs, it was decided to use mean as amplification factor and the correlation 
coefficient between RSA and PGA as weighing function. 

3.5.3 Spread in Predicted PGAs 

In order to study the spread in the predicted PGAs from SMAs and SRRs a 
parameter called average Coefficient of Variation was used. This is defined as the 
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average of coefficient of variation {aj/:) observed for all the stations across the 
predicted values and gives average extent of spread that occurs among the six generated 
PGAs at each station. Since the weighing functions will also be affecting the variation, 
the weight was assumed equal to 1.0 for this analysis. 

The coefficient of variation among the six calculated peak ground accelerations 
was computed for all the events (Table 3.9). It can be observed from the table that, 

(a) The spread of predictions in the case of SRRs is not very different from SMAs. The 
predictions based on SMAs differ from that based on SRRs by an average of 20% 
which ranges from 13% (Uttarkashi) to 50% (Kangra). 

(b) The spread in prediction is always more for SRRs compared to SMAs. 

It can be concluded that, the spread in PGAs calculated from SRRs and SMAs do not 
vary significantly. 

Based on these observations, it was decided to use the Mean values as conversion factors 
and CC as weighing functions for calculation of PGA from SRR recordings. The 
observed values of peak ground acceleration from SMAs for different events as well the 
predicted values of acceleration from the SRRs are given in Figure 3.12. It can be 
observed that the predicted values of acceleration seamlessly fit into the observed SMA 
data. 


• • • 
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Table 3.1 (a): Earthquake events considered in this studx 


Epicenter 


Lat‘N 


USGS 25.385 S 


25.200 S 


Location 

Array 

Date 

NE. India 

Shillong 

10Sepl986 

NE. India 

Shillong 

14 May 1987 _ 

NE. India 

Shillong 

6 Feb 1988 

NE. India 

Shillong 

6 Aug 1988 


9 Jan 1990 

N. India 

- 

:21 Aug 1988 
(Bihar-Nepal) 

N. India 

Uttar 

Pradesh 

20 Oct 1991 
(Uttarkashi) 

N. India 

Uttar 

Pradesh 

29 Mar 1999 
(Chamoli) 

N. India 

Kangra 

26 Apr 1986 
(Kangra) 


Source 


MD 


DE 


DE 


USGS 


DE 


use 5 


DE 


USGS 


USGS 


ISC 


USGS 


DE 


25.149 


25.392 


24.750 


26.750 


92.077 


91.600 



95.127 


94.533 


95.240 


86.620 




Table 3.1 (b): Earthquake events considered and number of recordings. 
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Table 3.2: Correlation of epicentral distances with average peak ground 
acceleration. 
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Table 3.3(a): Results of statistical analysis for the ratios of the data of 6 August 


Event 



iiSMA\ 

S3 1 

Mean 

SMA 

SRR 

2.157 

2.597 


4.903 

5.125 

2.109 

2.721 

2.016 

2.227 

4.458 

5.349 

1.338 

1.471 


IISRR 1 

104 1 

Median 

SMA 

SRR 

1.674 

2.438 

1.957 

1.856 

4.305 

4.496 

1.765 

2.438 

2.011 

1.842 

3.615 

4.761 

mtm\ 

1.337 

1.401 

1.252 

1.272 



2.157 2.597 1 1.674 2.438 1.635 1.276 0.713 0.777 - 1.415 


0.845 0.710 0.291 


0.631 0.611 - 0.303 


2.109 2.721 1.765 2.438 1.452 1.787 0.739 0.827 - 1.990 


2.016 2.227 12.011 1.842 0.710 1.759 0.874 0.716 - 0.994 


3.585 3.831 


0.121 0.525 


1.223 


• 2.391 



Table 3.3(b): Results of statistical analysis for the ratios of the data of 21 August 
1988 Bihar-Nepal earthquake. 
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Table 3.3(c): Results of statistical analysis for 
1991 Uttarkashi earthquake. 


the ratios of the data of 20 October, 


Event 


Uttarkashi 

20-Oct-91 



^MA 


Mean 


13 


#SRR 


Median 


95 


SMA 

SRR 

2.225 

2.422 

— Bl 

4.403 

3.587 

2.047 

2.580 

1.919 

1.578 

3.992 

3.508 

1.274 

1.518 

1.220 

1.556 

BEHMlHil 



1.890 1.375 



1.240 1.297 



CC 

SMA 

SRR 5 

0.820 

0.898 

0.943 

0.967 

0.701 

0.885 

0.879 

0.924 

0.947 

0.959 

0.747 

0.852 

0.991 

0.986 

0.998 

0.988 

0.997 

0.992 



Table 3.3(d): Results of statistical analysis for the ratios of the data of 29 March, 
1999 Chamoli earthquake. 
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Table 3.3(e): Results of statistical analysis for the ratios of the data of 26 April, 
1986 Kangra earthquake. 


Event 


^MA 


Mean 


2.238 2.654 1 1.972 2.824 0.911 1.405 0.820 0.972 -0.911 


2.606 1.808 2.583 1.533 1.027 0.865 0.859 0.946 1.988 


5.813 4.288 4.500 4.167 3.152 2.187 f 0.834 0.929 1.296 


\#SRR\ 

17 1 

Median 

SMA 

SRR 

1.972 

2.824 

2.583 

1.533 

4.500 

4.167 

1.970 

2.514 

2.460 

1.352 

4.290 

3.400 

1.218 

1.229 

1.182 

1.242 

1.115 

1.182 


5.049 1 7.139 I 4.290 1 3.400 | 2.140 10.070 0.903 0.603 -0.826 


. CC 

SMA 

SRR >5 

0.820 

0.972 

0.859 

0.946 

0.834 

0.929 

0.905 

0.891 

0.914 

0.767 

0.903 

0.603 

0.996 

0.885 

0.994 

0.978 

0.993 

0.871 
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Table 3.6: Spectrum amplification factors for horizontal elastic response (from 
Newmark and Hall, 1982), 


Damping 
% Critical 


.5 


One Sigma (84.1%) 



Table 3. 7: Correlation of amplification factors and epicentral distance for different 
earthquakes. 


Correlation factor Between Amplification Factor 
^■^residtnt I P^Aiverage > tind Epicentral Distance D 




- 0.418 


- 0.254 


- 0.105 


- 0.487 


- 0.252 


- 0.121 


20-Oct-1991 

Uttarkashi 

29-Mar-1999 

Chamoli 

+ 0.214 

+ 0.437 

- 0.358 

- 0.172 

- 0.343 

- 0.501 

+ 0.197 

+ 0.339 

- 0.437 

- 0.336 

- 0.349 

- 0.508 


- 0.373 


- 0.016 


- 0.327 


- 0.326 


- 0.107 


+ 0.528 


Table 3.8 Error (eqn 3.5) between PGAobsemdand PGApreduted calculated using 
different weighing factors and amplification factors (e.g. eqn 3.9). 
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Table 3.9 Coefficient of variation among the six PGA calculated from 5 !^ values of 
SMAs and SRRs (Weights equal to unity) 
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Figure 3.1: Tectonic features of the different regions considered in this study [modified 
from Rajendran et. al, 1992]. 




Figure 3.2: Relative position of arrays used in this study [modified fiom Chandiasekaran 
and Das, 1992]. 





Figure 3.3(a): Photograph of a SRR [Agarwal, 1991]. 
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T 


5” 0 Damping 


10 Damping 




Time Penod = 0 4 sec 



Sa = 9S (cro'sec') 


Sa= 8“ (cm,' sec') 



Time Period = 1 25 sec 



Sa = 3 b (cnv'sec*) 


Sa= 31 (cm'sec') 


Figure 3.3(b); Atypical SRR record [DEQ 1999]. 
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91.00 91.25 91.50 91.75 92.00 92.25 92.50 


Longitude (°E) 

Figure 3.5(a); Location of epicenters as reported by different agencies along with the 
recording stations activated for the Shillong earthquake of 10 September 
1986. 
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■ USGS Epicenter (M5.9) 

♦ DEQ Epicenter 

• SMA Station 


Longitude { E) 

Figure 3.5(b); Location of epicenters as reported by different agencies along with the 
recording stations activated for the Shillong earthquake of 14 May 1987. 
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■ USGS Epicenter (M5.8)’ 

♦ DEQ Epicenter 

• SMA Station 


Longitude ( E) 

igure 3.5(c): Location of epicenters as reported by different agencies along with the 
recording stations activated for the Shillong earthquake of 6 February 
1988 
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30.0 H 


■ USGS Epicenter (M7.3) 

♦ DEQ Epicenter 

• SMA Station 
O SRR Station 


28.0 - 


26.0 - 


24.0 - 


22.0 - 
Bl 


Longitude ( E) 


gure 3.5(d): Location of epicenters as reported by different agencies along with the 
recording stations activated for the Shillong earthquake of 6 August 1988. 
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92.0 92.5 93.0 93.5 94.0 94.5 95.0 95.5 


Longitude {°E) 

Figure 3.5(e): Location of epicenters as reported by different agencies along with the 
recording stations activated for the Shillong earthquake of 9 January 
1990. 
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75.5 76.5 77.5 78.5 79.5 80.5 

Longitude (°E) 

Figure 3.5(g); Location of epicenters as reported by different agencies along with the 
recording stations activated for the Uttar Pradesh (Uttarkashi) earthquake 
of 20 October 1991. 
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Tl.5 ■ 78.0 78.5 79.0 79.5 80.0 

Longitude (°E) 


Figure 3.5(h): Location of epicenters as reported by different agencies along with the 
recording stations activated for the Uttar Pradesh (Qiamoli) earthquake 
of 29 March 1999. 
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Figure 3.6: Schematic representation for the calculation of the resultant spectral 
acceleration from the SMA records. 
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Figure 3.7(e): Distribution of ratios of &4ri,f=io%/-^^r2,f=lO% October 20, 1990 

Uttarkashi event for SRRs and SMAs. 
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Ratio 
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IC=5%/C=10%, 0.75 



0 0.5 1 1.5 2 2.5 3 3.5 4 
Ratio 


0 0.5 1 1.5 2 2.5 3 3.5 4 
Ratio 


I (;=5%/^=10%, 1.25 


ir=5o/o/r=io%, 1.25 


0 0.5 1 1.5 2 2.5 3 3.5 4 


0 0.5 1 1.5 2 2.5 3 3.5 4 


Ratio Ratio 

(i) SMA (ii) SRR 

Figure 3.7(f): Distribution of ratios of (5^o.4scc,^=5%/‘^-^.4sec,f=lO% )> 

(5/40.75 sec,$-=5% / ‘5'-^0.75sec,$-=10%) (5'^115sec,$-=5%/‘S^1.25sec,g-=10%) 

for October 20, 1990 Uttarkashi event for SRRs and SMAs. 
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0 1 2 3 4 5 6 


Ratio 



Ratio 



0 1 2 3 4 5 6 

Ratio 


SRR 

Figure 3 .7(h): Distribution of ratios of SAj-^ ,g=WA / =10% August 21,1988 
Bihar event for SRRs. 
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Figure 3.7(1): Distribution of ratios of / ^^r 2 .f=l 0 % 26, 1986 

Kangra event for SRRs and SMAs. 
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0.5 1 1.5 

Rai 


(i) August 6, 1988 NE India earthquake 



1 1.5 

Ral 


(ii) October 20, 1991 Uttarkashi earthquake 
30 , 



0.5 1 1.5 2 2.5 3 

Ratio 


(iii) AugiJSt 21, 1988 Bihar earthquake, SRR 
Figure 3.8: Histograms for average values of ( 5 ^ 0 . 4360 , f=5%/'^).4sec, $-=10% )> 

(‘^'^.75seci^=5% / ‘^•^0.75 sec, 5'=10%) (‘^•^1.25sec,4-=5%/'^-^1.25sec,5-=10% 
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Figure 3.9: Schematic representation for the confidence interval for the null hypothesis, 
Ho. 
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Figure 3.10(a): Amplification factor of the resxiltant Spectral Acceleration versus 
epicentral distance for the August 6,1988 Shillong earthquake, best fit 
line used y = a\Q.x + b . 
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Figure 3.10(b): Amplification factor of the resultant spectral acceleration versus 
epicentral distance for the October, 1991 Uttarkashi earthquake, best fit 
line used y = a\n.x-¥b. 
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Figure 3.10(c): Amplification factor of the resultant spectral acceleration versus 
epicentral distance for the March, 1999 Chamoli earthquake, best fit line 
used y = alnx + b . 
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0 10 20 30 0 10 20 30 
Epicentral Distance (km) Epicentral Distance (km) 



Epicentral Distance (km) Epicentral Distance (km) 

Figure 3.10(d): Amplification factor of the resultant spectral acceleration versus 
epicentral distance for the April, 1986 Kangra earthquake, best fit line 
used y = alnx + b . 
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Figure 3.11: Variation of amplication factors for spectral acceleration from different 
attenuation relationships available in the literature =5% of critical). 
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Chapter 4 

Development of Attenuation Relationships 

4.0 Introduction 

This chapter studies attenuation of Peak Ground Acceleration (PGA) from nine 
Himalayan earthquakes using data from strong motion accelerographs (SMAs) as well as 
structural response recorders (SRRs). The observed PGA values are compared with the 
existing attenuation relationships available in the literature. Considerable differences are 
observed in the attenuation characteristics for different tectonic regions. New attenuation 
relationships have been developed based on Indian ground motion data. 

4.1 Indian Strong Motion data Versus Existing Attenuation Relationships 

The Seismo-tectonical environment of Himalayan region can be divided into 
three classes based on geology. They are (.1) The northern Himalayas (2) Central 
Himalayas and (3) North-Eastern Himalayas. Due to the constant thrust between Euasian 
and Indian plate (Chandra, 1992), the rocks in the northern Himalayas are in a fractured 
state. Gangetic foredeep in Central Himalayas is characterized by the presence of a deep 
alluvial plane. In North-Eastern Himalayas, the base rocks are considered to more 
homogenous compared to that in northern Himalayas. Presence of a subduction process 
has also been reported along the Burmese arc in this region. So, the earthquakes included 
in the study can be divided in to three categories based on geology of region of 
recording. (1) Earthquakes of northern Himalayas (1991 Uttarkashi M7.0 , 1999 Chamoli 
M6.6, and 1986 Dhaimshala Af5.5), (2) Earthquakes of Gangetic foredeep (1988 Bihar- 
Nepal M6.8), and (3) Earthqukes of North-East Himalayas (1986 NE M5.2, 1987 NE 
M5.9 , 1988 NE M5.8, 1988 NE Ml 3 , 1990 NE M6.1). Two events in the NE Himalayan 
earthquakes have their predicted epicenters near the Burmese arc. Accelerations 
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observed during the earthquakes are compared with the existing attenuation relationships 
to assess the suitability of their application in Indian region. 

The attenuation relationships used for comparison can be divided into three 
categories based on their origin. They are: 

(a) Relationships based on data from interplate earthquakes (Boore, 1997, Cambell, 1997 
and Sadigh et al., 1997) 

(b) Relationships based on data from subduction zone earthquakes (Young et al., 1997 
and Crouse, 1991) 

(c) Attenuation relationships based on data from intraplate earthquakes (Anders et al, 
1990; Toro et al., 1997). The intraplate earthquakes exhibit lesser rate of attenuation 
compared to most of the attenuation relationships which are mostly developed based 
on the data from western America, where the base rocks exist in a highly fractured 
state. - 

Since different relationships require different types of parameters as input, the 
following assumptions were made for the parameters used in them for Indian 
earthquakes. 

(1) . The site conditions are assumed rocky (Shear wave velocity = 620m/s ) for 
all events except for Bihar-Nepal event. Since the Bihar-Nepal earthquake originated in 
the Gangetic alluvium (Shear wave velocity = 310 m/s) propagation medium was soil. 

(2) The information about the source mechanism of the most of the Indian 
earthquakes are not available. Hence, the source mechanism was assumed be reverse 
thrust faulting, for all events expect for the 1988 Bihar-Nepal M6.8 earthquake. It was 
observed that assumption of source mechanism as reverse fault predicted highest 
accelerations compared to all other mechanisms. The 1988 Bihar-Nepal M6.Z earthquake 
occurred in deep alluvium and fault mechanism identified was of strike-slip nature. 
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(3) Since the information about length and direction of fault was not available, 

the closest distance to^ surface projection of fault was assumed equal to epicentral 
distance. 

(5) The closest distance to plane of mpture was assumed to be equal to 
hypocentral distance. 

(6) Youngs et al, (1997) have used the geometric mean of PGA as the predicted 
parameter. The maximum value of two horizontal components of PGA was used by 
Anders et ah, (1990) as the output parameter. But, these relationships have not provided 
any comparison of these parameters with respect to the average PGA, which is the 
parameter used in this study. So wherever these equations were used for predicting the 
PGA, the result was used without any modification. 

The observed data points from different events in the Himalayan region and the 
predicted acceleration are shown in Figures 4.1 to 4.3. For ease of comparison the plots 
have been made based on their origin. The least square fit of each earthquake is also 
plotted to compare the deviation of the predicted data from the ideal relationship. Based 
on the plots, the following general observation can be made: 

1. Interplate relationships give a better representation of the observed data for 
earthquakes in north & North-Eastern Himalayas compared to intra plate 
relationships. 

2. The accelerations predicted by the interplate relationships of Campbell, (1997) and 
Sadigh et al., (1997) underestimate the accelerations observed from the earthquakes 
in the North-Eastern Himalaya. The difference in predictions is maximum for the 
earthquakes generated from Burmese arc. Of all the relationships used for 
comparison, Boore (1997) seems to give a better representation of the data. 
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3. Two North East events namely 1988 NE Ml. 3 and 1990 NE M6.1, whose epicenters 


were situated along the Burmese border, generated a much higher acceleration 
compared to the other events reordered in the region. For these events (1988 NE 
Af7.3and 1990 NE M6.1), the predictions by Youngs et al, (1997) are found to be 
closer to the observed values compared to the relationship developed by Crouse, 
(1991). For M7.3 event, the rate of attenuation of PGA is found to be less than that 
predicted by these relationhips. 

4. For North India Earthquakes namely 1986 Dhaimshala M5.5, 1991 Uttarkashi 3/7.1 
and 1999 Chamoli 3/6.6, the existing interplate relationships are able to predict the 
accelerations without much error. But the optimum relationships are found to be 
different ones for each event. For Chamoli and Kangra events, the intraplate 
relationships predicted much higher accelerations compared to the observed values. 

5. The rate of attenuation predicted by the intraplate relationships is in line with the 
attenuation observed in 1988 Bihar-Nepal 3/6.8 earthquake. The rate of attenuation 
of PGA for 1988 Bihar-Nepal 3/6.8 appears to be very rapid. An extrapolation of the 
observed acceleration indicates that the peak ground acceleration may have exceeded 
Ig near the source of the earthquake. 

6. The North-East event with magnitude 6. 1 has produced very high accelerations and 
all the relationships under-predict the accelerations. 

The deviation of the points from the observed data was calculated by averaging 
the squares of their individual deviations from the observed values in logarithmic scale. 
The square root of the average can be used as a measure of average distance from the 
predicted curve and the observed data, i.e.. 


Error = 




) ln(y predicted ))^ 


j=l 


(4.1) 
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The results are given in Table 4.1. Most of the deviations observed in the data are greater 

than the standard deviation of the models considered. For the 1987 NE M5.2 1988 NE 

and M5.9 events, the calculated deviation 'from the prediction was found to be close to 

the standard deviation predicted by the attenuation relationship of Boore et al. (1997). 

For all other events, the sample points were generally found to lie more than ±la-away 

from the predicted points. The relationship developed by Sadigh et al. (1997) and 

Campbell, (1997) gives a very high average error of prediction when compared with all 

the observed events. Thus, the earthquakes under consideration cannot be assumed to 

come from the same parent population as that used for the development of the various 

attenuation relationships. 

The comparison of the Indian strong motion data with some of the existing 
attenuation relationships shows that in many cases the existing attenuation relationships 
have failed to represent the observed PGAs and a new attenuation relationship has to be 
generated for this region. The new attenuation relationship, which will be based on the 
enlarged database generated through the current work, can be used for prediction of the 
acceleration for future Himalayan earthquakes. 

4.2 Development of Attenuation Model 

4.2.1 Selection of Parameters 

The parameter, which is being predicted, is known as dependent variable and the 
parameters used for the prediction are called independent variables. In the present study, 
average PGA is taken as the dependant variable. 

4.2.1. 1 Earthquake Parameters 

Seismic moment is one of the ideal parameter, which can be used to represent the 
earthquake size. Because of the lack of information on the seismic moment of the 
earthquakes under consideration, earthquake magnitudes mi, and Ms are used in the 
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present study for characterizing the size of an earthquake. Magnitudes used in the 
analysis can be considered to be equivalent to moment magnitude because, below 

THb is considered to be equal to My, and above jnb~6.0, Ms is considered to be 
equal to the moment magnitude. 

4.2.1.2 Source Parameters 

In the present study, all the events occurred near the plate boundary and are 
interplate earthquakes. Since the fault mechanism of all the events are not well defined, 
the style of faulting is not included as a variable in the study. 

4.2.13 Propagation Parameters 

The propagation parameters characterize the effects due to wave scattering, 
inelastic attenuation and geometric attenuation of acceleration. Since the type and length 
of faulting, amount of slip etc. are not available, the distance measurements such as 
distance to horizontal projection of fault mpture, closest distance to mpture, distance to 
zone of maximum energy release etc. could not be used as the propagation parameter. In 
this analysis epicentral distance is used to represent the propagation parameter. 

4.2.1.4 Site Parameters 

Other than the broad classification on the local site conditions [Chandrasekaran 
and Das, 1994], no detailed information is available about the local site geology. Due to 
the lack of this information, all the events considered in this study, except the earthquake 
in Gangetic alluvium (1988 Bihar-Nepal M6.8), are assumed to be recorded on rocky 
strata since it is the general characteristic of the region. 

4.2. 1.5 Data Selection 

The error in the predicted values of a parameter depends on the type of data that 
is used in the regression analysis. The scatter in the data to be analyzed is minimal if the 
recordings are taken from a tectonic region of similar attenuation-source features. 
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consistent record processing techniques are employed and recording instruments of 

similar dynamic characteristics [Campbell, 1985] are used. The database was assumed to 

be a homogenous one and no initial segregation of the data is done. 

A plot of the epicentral distance V 5 . magnitude of the event was made for all the 
observations (Figure 4.4). It was found that a Magnitude of the earthquake and the 
distance to the farthest station follows a particular relationship i.e., there is a cutoff 
distance associated with each magnitude. For the data set used in the study, the cut-off 
line can be approximately denoted by M= 4.8 + 0.0038J?, where R is the epicentral 
distance. 

4.2.2 Selection of Attenuation Model 

After deciding the data set as well as the group of independent variables that will 
be used in the analysis, it is necessary to select a proper attenuation model for the data. If 
the model is to be used for prediction using a data, which lies inside the range of 
parameters used for regression analysis, then a simple attenuation model can be used. In 
order to extrapolate the data from the observed ones, the model should have a physical 
basis during formulation [Campbell, 1985]. 

The amplitude of motion of a harmonic wave in an infinite elastic half space is 
given by [Dahle et al., 1990] 

(4.2) 

where A is the observed ground motion amplitude, R is the hypocentral distance, m is 
the magnitude and Aq, a, b and q are constants. The linearisation of Eq. (4.1) leads to 

In.^ = In^ +am -hlni? + ^i? . (4.3) 

Considering this as the basic form of attenuation, more variables can be added to 
represent the dependence between magnitude and distance, site parameters etc. So, the 
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attenuation relationship that can be used in a region can be of the form [Campbell, 1985] 

Y = bi.fi(m)f 2 {R)fiim,R)f^P, ^ (44) 

where Y is the strong motion parameter predicted; usually takes the form 

fl (^) = • This is based on the definition of magnitude as a logarithmic variation of 

amplitude [Richter, 1956]. Sometimes the quadratic form of magnitude is also used 
[Toro et al, 1997; Spudich et al, 1997], where the expression takes the form 

f\ i^) = ci{m-k) + C 2 (m-kf (4.5) 

where A: is a constant. During analysis, k is selected in such a way that it gives the 
minimum squared error. A cubic term of magnitude has also been used, which is as 
follows [Joyner and Boore, 1991], 

/l (m) = Ci(m-k) + C 2 (fn-k)^ + C 3 (m - kf . (4.6) 

A third order equation was used instead of linear or second order equation, after it was 
found that the use of cubic term resulted in a better fit. The commonly used form of 
distance related term /jCT?) is given by 

f2(R) = e'’*‘‘(R^■bsr'^. (4.7) 

The exponential term used in this equation represents the material damping and 
scattering, which is also known as inelastic attenuation of ground motion. The second 
term denotes for geometrical spreading of energy. The geometrical attenuation rate is 
given by the coefficient 63 [Joyner and Boore, 1991]. The coefficient 65 is used to limit 

the value of Y at zero distance. This becomes more significant when the distance 
parameter being used is epicentral distance or shortest distance to the surface projection 
of rupture. Another widely used form of / 2 (R) [Ambraseys, 1995; Joyner and Boore, 
1981; Fukushima and Tanaka, 1990] is given by 
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The constant term in the brackets 'is considered analogous to hypocentral distance. In 
order to account for the combined effect of magnitude and distance, the form of the 
function used is given by [Campbell, 1997; Sadigh et al, 1997] 

/3(m,^) = (R-66£'^"')■'’. (4,9) 

Depending on the effect of local site conditions, the style of faulting, soil structure 
interaction parameters and type of earthquake the variables, which can take pre-assigned 
values, can also be introduced in the formulation. This is done through the introduction 
of dummy variables corresponding to these parameters and is denoted through the term 
UP. 

4.3 Attenuation model for Himalayan Earthquakes 

The number of well-recorded earthquakes available for the Himalayan region is 
limited. So, the use of an exhatistive model including all the parameters mentioned above 
should be applied with caution. A simple form of attenuation equation compared to the 
one discussed above was chosen because of (a) lack of enough near source recordings, 
and (b) inadequate number of earthquakes to arrive at a meaningful nonlinear 
relationship with respect to magnitude, z.e., the model contains only linear term in 
magnitude. In R and R terms to incorporate the effect of attenuation due to geometric 
spreading as well as material damping, and a soil dependant parameter. The form of the 
attenuation model is ^ven as 

\xiY = b\ +^>2-^ + b2,R'^b^\xiR-^b^H + {b(^ ■¥bjR)S . (4.10) 

The term if is added to include the effects of higher acceleration (if any) observed 
in the subduction zone earthquakes. The value of /f is 1.0 for the earthquakes subduction 
zone earthquakes (1988 NE Ml. 5 and 1990 NE M6.1), and 0.0 for all other earthquakes. 
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A dummy value S is introduced to consider the change in the rate of attenuation in 

alluvial medium. S is equal to 1.0 for 1988 Bihar-Nepal A/6.8 earthquake and for all 
other events which had rock as medium of propagation, S is zero. 

4.3.1 Preliminary analysis of data 

Before the data is regressed to get a particular attenuation model, it is necessary 
to do some preliminary analysis on the data to eliminate inconsistent data as well as to 
find the form of variables that gives a better fit for the data. The plot of PGA V 5 
epicentral distance of the data used in the study is given in Figure 4.5. 

The regression analysis for the model given by Eq. (4.10) was carried out using 
MATLAB subroutine REGRESS and G02DAF fi-om NAG library. Both the subroutines 
use method of least squares for the regression analysis. The results of regression using 
the entire data of nine earthquakes are shown in Table 4.2. The table also gives the 
standard error, which denotes the amount of uncertainty in the calculated coefficients. 
This is generally not reported along with the attenuation relationships, but is useful for 
understanding the level of confidence that can be put on the predicted coefficients. T hi s 
part can be dealt in detail after the regression form has been finalized. 

The coefficient of earthquake magnitude is quite low indicating a very low 
dependence on magnitude or due to some error related to the sample considered in the 
analysis. For the model under consideration, a change of magnitude from 5.0 to 6.0 
increases acceleration only by 12%, which is very less compared to the difference in 
amount energy released in the earthquake. For an increase in magnitude by 1.0, the 
average increase in acceleration predicted in the literature [e.g., Boore et ah, 1997; 
Spudich et ah, 1997; Ambraseys, 1995; Fukushimaand Tanaka, 1990] is about 50% and 
it ranges fi-om 25 to 75%. This anomaly may be attributed to the inconsistency of 
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observations in the data used, inability of the model to represent the complete database or 
the true characteristic of the region. 

4.3.2 Study of Outliers 

Any experimental data recorded will contain some outliers, which are to be 
removed before the data is analyzed. These can be either erroneous measurements, which 
are not consistent with the other observed points or points that get undue importance 
during analysis. When there are two or fewer explanatory variables, scatter plots will 
quickly reveal any outliers in X. When number of coefficients are greater than 2 , scatter 
plots may not reveal multivariate outliers, which are separated from the bulk of Z-points 
but do not appear as outliers in a plot of any single explanatory variable or pair of them. 
A least square projection matrix known as hat matrix [H] [Belsley et al, 1980] can be 
used as an indicator of the influence of the observation on the predicted value. The hat 
matrix is defined as 

H = x)' X'^ , (4.11) 

where X is the nxp matrix containing explanatory variables. Hat matrix can be used to 
calculate the fitted predicted values, since 

y = Xb = Hy. (4.12) 

The contribution of each observation y,- in the prediction of y/ is obtained from the hat 
matrix. The influence of the response value, y,- on the fit is most directly reflected in its 
impact on the corresponding fitted value y,- . This information can be obtained from the 

diagonal elements of H, known as {A} vector. This {h) vector can be used as a tool to 
obtain the information about the over influencing observations. The proof for the fact 
that the average size of the diagonal element A,- is pjn and Ipjn provides a good cut-off 
for the removal of the influential data is available in literature [Belsley et al, 1980]. 
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For the model given in Eq. (4.9), the values of {h} vector as computed are plotted 

with respect to epicentral distance of the observed points (Figure 4.6). It is seen that for 
each earthquake, the plot of {h} vector follows almost a parabolic shape. It is evident 
from the plots that the near source recordings (<10km) by SMAs as well SRRs are highly 
influencing the predicted values. The recordings at distances more than SOOkms were 
also found to be highly influential in prediction of the data and these were records 
corresponding to August 6, 1988 earthquake only. Majority of the points from 1988 
Bihar-Nepal M6.8 earthquake is found to fail in satisfying the stipulated criteria. The 
failure of large number of points corresponding to a particular earthquake indicated a 
different behavior being exhibited by these points compared to the other events. The 
1988 Bihar-Nepal M6.8 earthquake was the only event that was recorded after 
propagating through a deep alluvial plane. All other events had rock as the medium of 
propagation. Subsequently, 1988 Bihar-Nepal M6.8 earthquake data was removed from 
the database to be separately regressed later. The new database was regressed with the 
model without soil dependant coefficients as given below, 

]iiY = bi+b2M+b2R + b4\nR + b5H . (4.13) 

The results are given in Table 4.3. The magnitude dependent coefficient did not 
show any marked improvement even after removal of Bihar-Nepal earthquake. This 
model also indicates that there are still data points, which have strong influence on the 
model. 

The analysis of {h} matrix indicates that, the influential observations are detected 
only at either extreme (with respect to epicentral distance) of the observations. One 
influential row from the matrix containing explanatory variables was deleted at a time 
and the values of {A} were re-computed. In general the influential observations consisted 
of recordings at distances less than 10km and recordings at distances more than 500 km 
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(1 recording from 1986 NE M52, 14 from 1988 NE AO.S, 12 from Uttarkashi, 6 from 
Chamoli and 4 from Kangra events). The regression model after removal of the 
influential observations, showed an even weaker dependence (decreased to 8 % from 
12 %) of peak ground acceleration on magnitude (Table 4 . 3 ). The coefficient 65 , which 

was introduced to incorporate the difference in attenuation characteristics shown by the 
earthquakes from Burmese arc (subduction zone), was found to have a positive sign 
indicating magnification of the accelerations (about 300%) compared to a similar 
magnitude earthquake from a different source. Since the problem of magnitude 
independence persisted even after the removal of the influential data and the regression 
form did not seem to exclude any important parameter, the data was checked for any 
difference in characteristics, i.e., the low magnitude dependence of the equation might be 
due inclusion of data with different characteristics in the same regression form. 

The plot of PGA versus epicentral distance for all the earthquakes after removal 
of influential observations are shown in Figure 4.7. Figure 4.8 shows the PGA recorded 
for different events in the North East India and Figure 4.9 shows those for the North 
India events. Figure 4.10 shows the variation of PGA for Bihar-Nepal earthquake. 
Individual least-square fit lines of all the events are plotted and are given in Figure 4.1 1. 

4.33 Observations on variation of PGA of different events: 

Earthquakes recorded in the North-East India seem to have similar acceleration 
when compared to events with higher magnitude in North India. From Figure 4.1 1 it can 
be seen that an earthquake with a magnitude 5.2 in the North-East Himalayas has 
produced comparable acceleration with an earthquake of magnitude 6.6 (Chamoli) in 
northern Himalayas. Similarly, an earthquake with magnitude 5.8 in North-East 
Himalayas has produced similar acceleration for an earthquake of magnitude 7.0 
(Uttarkashi) in northern Himalayas. Also, the acceleration recorded during Kangra 
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earthquake is lower than that of magnitude 5.2 earthquake in North-East. If a regression 
analysis is conducted without considering these characteristics, because of apparent 
independent behavior of acceleration with respect to magnitude, the magnitude related 
regression coefficient tends to become very low. So in order to get a meaningful 
magnitude dependant relationship, the separation of the database based on the region is 
required. Since the subduction zone events were found to produce higher acceleration, it 
was decided to have separate regression form for these events. 

Though there is a large difference in earthquake magnitudes, the PGA generated 
by 1990 NE M6.1 are comparable with that produced byl988 NE Ml 3 earthquake. This 
anomaly can be attributed either to an error in magnitude measurement or error in the 
epicenter calculation. The average error in the calculation of magnitude by different 
agencies is ±0.4 (Table 3.1a). Even if this variation is taken into consideration, there 
will be a difference of 0.8 in magnitude and the similarity in the order of accelerations 
cannot be explained. 

It is also evident from Figures 3.5(a) to 3.5(d) diat there is an error of 50-100km 
in locating the epicenters of North-East India earthquakes by different agencies. For 
events 1986 NE M5.2, 1987 NE.V/5.9, 1988 NE M5.8, and 1988 NE M7.3, the epicenters 
calculated using the strong motion data were found to lie closer to the recording site 
compared to the epicenter predicted by USGS. For January 20, 1990 earthquake, the only 
available location of the epicenter is from USGS. Thus, the best epicenter cannot be 
chosen from different locations available based on the correlation of distance and 
acceleration. There is a possibility that the error in locating the epicenter may have crept 
into the 1990 NE M6.1 event Under these circumstances, it becomes necessary to 
remove the 1990 NE M6.1 event from the database. 
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Based on the above observations, the database can be classified into 4 categories 
corresponding to four regions. They are, 


(a) North India earthquakes. 


(b) North East India non-subduction zone earthquakes, 

(c) NE India subduction zone earthquakes, and 


(d) Earthquakes generated in the Gangetic foredeep, which is basically situated on 
alluvial plane. 


4.4 Development of Attenuation Relationships 

Tables 4.4 to 4.8 gives the coefficients of regression related to each region along 
with a of the model. It can be observed from Table 4.4 that the regression predicts a 
negative value of material damping for North East India non-subduction earthquakes. 
Thus, the regression analysis for this region is carried out after removal of material 


damping term (b^R) from the model. The resulting regression coefficients along with 


the Standard Error (SE) in prediction of coefficients and SD of the model are given in 
Table 4.5. It can be observed that, the magnitude term has come out to be a significant 
value when the regression analysis was done separately for each region (Table 4.5 to 
4.6). 

For subduction zone event and earthquakes in alluvial planes, only single event is 
available. This necessitated the regression analysis of the data without magnitude term. 
The model used for regression analysis is given as 

\nY = bi+b2R + b3]nR (4.14) 

The results of regression analysis of these two groups are given in Tables 4.7 and 
4.8. For the 1988 Bihar-Nepal M6.8 earthquake, both the coefficients corresponding to 
material damping as well as geometric spreading terms are higher compared to August 6, 



If7.0. Eq.4.17 for the subduction zone earthquake is for M7.3 and Eq.4.18 for the Bihar- 
Nepal region for M6.8 event. To plot equations for an Af7.0 event for these two 
categories, the b 2 value was assumed as 0.70 and the value of bi adjusted accordingly. 

Considering the epicentral distance range of 10km to 100km, which is of 
maYimiim engineering interest, motion is weakest in the Central Himalayan earthquakes 
and strongest in the Bihar-Nepal earthquake. For the same magnitude {MIS)), the ratio of 
PGA in the Indo-Gangetic plains to that in the Central Himalaya ranges from 5.2 (at 
10km) to 3.4 (at 100km). The rate of attenuation with distance is lowest for the 
subduction event in North-East India and is highest for the Bihar-Nepal earthquake; ratio 
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of PGA at 10km to that at 100km is 2.6 for the former and 9.9 for the latter. Since the 

number of near source recordings was less, they^were not included in the regression 
analysis. So this attenuation relationships may not be applicable in the near source 
region. 



The adequacy of the model can checked either by analysing the residuals or by 
testing the significance of each coefficient derived during the analysis. 


4.5.1 Analysis of Residuals 

The overall adequacy of the model is best assessed from an analysis of residuals. 
By the very nature of regression analysis, the residuals have a mean close to zero. If no 
apparent trend in the residuals is observed, then the model can be considered adequate 
[Draper and Smith, 1998]. A trend indicates an inadequacy in the model to predict the 
data and would require modification of the functional form. The residuals plotted with 
respect to epicentral distances (R) are given in Figures 4.14 to 4.17. No particular trend is 
shown by any of the plots from which we can ascertain that the models used in the 
regression analysis are adequate to represent the variation in data. 

4.5.2 Significance Testing of Coefficients 

The significance of the calculated coefficients can be estimated by means of a t- 
test. The t statistic is defined as [Draper and Smith, 1998]: 


t = 


SEih) 


(4.19) 


Under the normality assumption, the variable bi follows t distribution with N-P degree 
of freedom, where P is the number of coefficients in the equation. If the true value of 
is specified in the null hypothesis, the t statistic can be computed from the available 
sample and used as a test statistic. 
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All the coefficients calculated are tested with a level of significance of 5% for the 

particular value being equal to zero. The null hypothesis being Hq : = 0 against the 

alternate hypothesis Hi'. 0 . The calculated t- statistic as well as t^ji n-P for the 

data used in the analysis is given in Tables 4.5 to 4.8. All coefficients, which are 
calculated for North India as well as North East India non-subduction zone earthquakes, 
differ significantly from the assumed 0 with 5% level of significance. For North East 
Indian subduction zone earthquake (1988 NE Af7.3) as well as for the 1988 Bihar-Nepal 
Af6.8 earthquake, the null hypothesis of these coeflScients being equal to zero cannot be 
rejected when tested with a level of significance of 5%. This implies, that the values 
obtained for these events should be used with caution. 


• • • 
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Table 4.1: Average error of prediction from the observed data for different attenuation 
relationshivs. 


Events \ Boore \Sadigh\Campbell\Youne\Crouse\ Ander,\ Tnm 


lO-Sep-86 0.5425 1.6036 0.9887 0.5043 0.6635 0.6864 |0.7157 


1 8-May-87 0.6132 2.4032 1.6330 1,0438 1.0846 " 1.3963 j 12375 

06-Feb-88 0.4712 1.9119 1.2179 0.7449 0.8273 0.9456 iQ.8695 

06-Aug-88 0.8678 2.4485 1.8430 1.0527 0.7943 2.2368 j 1.571 Q 

20-Jan-90 1.2823 3.7612 2.7534 2'.0318 1.8440 3.0760 i2.3614 

20- 0ct-90 0.7160 1.0881 0.6726 0.6276 0.6438 0.7547 j0_5923 

29-Mar-99 0.7856 1.0988 0.8915 0.9294 0.9014 0.9148 il.0105 

26-Apr-86 0.6993 0.8122 0.6924 0.9745 0.8331 0.8018 lQ.8138 

21- Aug-88 0.9957 2.2212 1.7076 1.0168 0.9071 1.5284 il.3054 

Table 4.2: Coefficients of regression for the model given below using the database 
containing accelerations from all the events. 

+ &3 i? + In i? + b^H + {b^ + b^R^S 

Value Mean Standard 
Error 
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Table 4.3: Coefficients of regression for the model given below using the database 
after the removal of data obtained from 1988 Bihar-Nepal M6.8 
earthquake. 

lnY = bj+b2M +b3R + b4hR+b5H 


■ 

Data with all points 

Data after removal of 
outliers 

Value SE 

Value 

SE 

h 

-1.545 

0.4188 

-2.559 

0.5751 

h 

0.1176 

0.0657 

0.0629 

0.068 

h 

- 0.002 

0.0005 

-0.005 

0.0009 

h 

-0.511 

0.0707 

-0.145 

0.136 

h 

1.2429 

0.0804 

1.2986 

0.0814 

SD 

0.5887 

0.5855 


Table 4.4: Coefficients of regression for the model given below using the accelerations 
recorded by 1986 NE M5.2, 1987 NE M5.9 and 1988 NE M5.8. 

In y = bi b 2 M + b^R + In J? 


Coef. 

Regression values 

Value 

SE 

h 

-3.1373 

1.8548 

bi 

0.7265 

0.3362 

h ' 

0.0013 

0.0050 

h 

-0.9515 

0.4950 

SD 

0.4478 1 


Table 4 5" Coefficients of regression for the model given below using the database 
accelerationsrecordedbyl986NEM5.2, 1987 NE 5.9 and 1988 NE MS. 
\xiY = bi+b2M-t-b3\nR 

Niimher of observations = 43, 40 =2.02. 


of 

8 . 


Coef. 

Regression values 

Value 

SE 

t Statistic 

h 

-3.441 

1.4381 

-2.393 

h 

0.706 

0.3236 

2.182 

h' 

-0.828 

0.1625 

-5.096 

SD 

0.4368 
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Table 4.6: Coefficients of regression for the model given below using the database of 
accelerations recorded by Uttarkashi, Chamoli and Kangra earthquakes. 

In F + b 2 M + b^R + 6^ In 


Number of observations =139, ^ 50 /^^ 135 = 1.97. 


Coef. 

Regression values 

Value 

SE T Statistic 

h 



-4.888 

h 

0.5856 

0.1524 

3.843 

h 

-0.0032 

0.0019 

mm 

■■ 



-2.009 

SD 

0.5973 


Table 4. 7: Coefficients of regression for the model given below with only subduction 
zone event (1988 NE M7.3). 

\nY = bj+ b 2 R + b^ Ini? 

Number of observations = 123, ^ 50 / 0 ^ 120 “ 


Coef. 

Regression values 

Value 

SE 

t Statistic 

h 

-0.420 

1.5987 

-0.262 

h 

-0.004 

0.0017 

-2.294 

h 

-0.241 

0.3677 

-0.655 

SD 

0.5776 


Table 4.8: Coefficients of regression for the model given below with only 1988 Bihar- 
Nepal M6.8. 

]xiY = bi+b2R + bs\nR 


Number of observations = 38, ~ 


Coef. 

Regression values 

Value 

SE 

t Statistic 

h 

2.1034 

3.1041 

0.6776 

bl 

-0.006 

O.OOM 

- 1.111 

63 

-0.76 

0.7881 

-0.965 

SD 

0.6968 1 
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Table 4.9: Average error of prediction from the observed data for different attenuation 
relationships. 


Events Boore Sadigh Campbell Young\ Crouse \Anders 


10-Sep-86| 0.541 1.662 1.032 


18-May-87 0.613 2.403 1.633 


06-Feb-88 0.471 1.912 1.218 


06-Aug-88 0.895 2.389 1.825 


20-0ct-90 0.721 1.094 0.676 


29-Mar-99 0.824 1.127 0.923 


26-Apr-86 0.827 0.941 0.799 


21-Aug-88 0.996 2.221 1.708 


E! 


S 

Ifes 


0.717 


1.396 


0.946 


0.761 


Toro 


0.695 


1.238 


0.870 


86 


0.596 


New 


0.435 


0.372 


0.449 


0.959 0.907 


0.958 

1.034 

0.715 

0.939 

0.911 

0.633 

1.528 

1.305 

0.669 
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Figure 4.1(a); Distribution of peak ground acceleration for 10 September 1986, NE India 
earthquake and the predicted values of PGA using different relationships 
available for interplate earthquakes. 



Figure 4.1(b): Distribution of peak ground acceleration for 14 May 1987, NE India 
earthquake and the predicted values of PGA using different relationships 
available for interplate earthquakes. 
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Figure 4.1(c): Distribution of peak ground acceleration for 6 February 1988, NE India 
earthquake and the predicted values of PGA using different relationships 
available for interplate earthquakes. 



Figure 4.1(d): Distribution of peak ground acceleration for 6 August 1988, NE India 
earthquake and the predicted values of PGA using different relationships 
available for interplate earthquakes. 
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Figure 4.1(e): Distribution of peak grotmd acceleration for 10 January 1990, NE India 
earthquake and the predicted values of PGA using different relationships 
available for interplate earthquakes. 



Figure 4.1(f); Distribution of peak ground acceleration for 20 October 1991, Uttarkashi 
earthquake and the predicted values of PGA using different relationships 
available for interplate earthquakes. 
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Figure 4.1(g): Distribution of peak ground acceleration for 29 March 1999, Chamoli 
earthquake and the predicted values of PGA using different relationships 
available for interplate earthquakes. 



Figure 4.1(h): Distribution of peak ground acceleration for 26 April 1986 Kangra 
earthquake and the predicted values of PGA using different relationships 
available for interplate earthquakes. 
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Figure 4.1(i): Distribution of peak ground acceleration for 21 August 1988, NE India 
earthquake and the predicted values of PGA using different relationships 
available for interplate earthquakes. 




1 10 100 1000 
Epicentral distance (km) 

Figure 4.2(a): Distribution of peak ground acceleration for 10 September 1986, NE India 
earthquake and the predicted values of PGA using different relationships 
available for intraplate earthquakes. 



Epicentral distance (km) 

Figure 4.2(b): Distribution of peak ground acceleration for 14 May 1987, NE India 
earthquake and the predicted values of PGA using different relationships 
available for intraplate earthquakes. 
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Figure 4.2(c): Distribution of peak ground acceleration for 6 February 1988, NE India 
earthquake and the predicted values of PGA using different relationships 
available for intraplate earthquakes. 



Figure 4.2(d): Distribution of peak ground acceleration for 6 August 1988, NE India 
earthquake and the predicted values of PGA using different relationships 
available for intraplate earthquakes. 
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Figure 4.2(e): Distribution of peak ground acceleration for 10 January 1990, NE India 
earthquake and the predicted values of PGA using different relationships 
available for intraplate earthquakes. 



Figure 4.2(f): Distribution of peak ground acceleration for 20 October 1991, Uttarkashi 
earthquake and the predicted values of PGA using different relationships 
available for intraplate earthquakes. 
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Figure 4.2(g); Distribution of peak ground acceleration for 29 March 1999, Chamoli 
earthquake and the predicted values of PGA using different relationships 
available for intraplate earthquakes. 



Figure 4.2(h): Distribution of peak ground acceleration for 26 April 1986 Kangra 
earthquake and the predicted values of PGA using different relationships 
available for intraplate earthquakes. 
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Figure 4.2(i); Distribution of peak ground acceleration for 21 August 1988, NE India 
earthquake and the predicted values of PGA using different relationships 
available for intraplate earthquakes. 
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Epicentral distance (km) 

Figure 4.3(a); Distribution of peak ground acceleration for 6 August 1988, NE India 
earthquake and the predicted values of PGA using different relationships 
available for subduction zone earthquakes. 



Epicentral distance (km) 

Figure 4.3(b): Distribution of peak ground acceleration for 10 January 1990, NE India 
earthquake and the predicted values of PGA using different relationships 
available for subduction zone earthquakes. 
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Figure 4.4; Scatter-graph of the epicentral distances from stations and magnitude of each 
events recorded by them. 
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Figure 4.5; Plot of the variation of accelerations recorded in the different earthquakes. 

with respect to the epicentral distance. This includes database consisting 
accelerations from all events. 




Values of hat matrix, hi 
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Figure 4.6: Distribution of k with respect to epicentral distance for the data set 
considered in the analysis. The dashed line gives theoretical cutoff for the 
data. 
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Figure 4.7: Plot of the variation of accelerations recorded in the different earthquakes. 

with respect to the epicentral distance. This includes database consisting 
accelerations from all events after removal of outliers. 
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Epicentral Distance (km) 


Figure 4.8: Plot of the variation of accelerations recorded in the NE Indian earthquakes 
with respect to the epicentral distance. 
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Figure 4.9: Plot of the variation of accelerations recorded in the N Indian earthquakes 
(Uttarkashi, Chamoli and Kangra). 
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Epicentral Distance (km) 

Figure 4.10: Plot of the variation of accelerations recorded in the Bihar- Nepal 
earthquake with respect to the epicentral distance. 
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Figure 4.11; Plot of the best fit line of the accelerations recorded for each events with 
respect to the epicentral distance. 
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Proposed, Boore, Sadiqh, Campbell, Crouse, Younas, SMA,^ SRR. 

Figure 4.12; The predicted values of acceleration using the new attenuation relationships 
plotted for different events along with the existing attenuation relationships. 



PGA 'g' 
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Figure 4.13 Comparison of the proposed equations for the four tectonic categories for a 
M7.0 event 




Epicentral distance (km) 

Figure 4.15: Plot of residues of the regression analysis for NE, subduction zone 
earthquake (August 6, 1988). 






Figure 4.16; Plot of residues of the regression analysis for Central Himalayan 
earthquakes. 
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Figure 4.17; Plot of residues of the regression analysis for Bihar-Nepal earthquake 



Chapter 5 


Summary and Conclusions • 


From pre-historic times, earthquakes are considered as one of the most dreaded 
natural hazards. Apart from loss of human lives, earthquakes cause a lot of damage to the 
region in terms of collapse of dwellings, lifelines and loss of producti\ity. However, it is 
one of the few hazards, where humanity is winning the race in terms of the minimization 
of damages caused. This is achieved through the proper methodology followed in 
analysis, design, detailing and construction of structures. Estimation of reasonable level 
of peak ground acceleration for which the structure should be designed during an 
earthquake is an important step in the design process. 

The available information about the attenuation of peak ground acceleration in 
Indian earthquakes is very limited. During the recent earthquakes in India, a large number 
of recordings from Structural Response Recorders have been generated as well as 
recordings from Strong Motion Accelerographs. The database of peak ground 
acceleration on Indian earthquakes can be increased if a proper method can be identified 
to relate the responses from structural response recorders to peak ground acceleration. 
The aims of the present study are as follows: 

- Study the responses from Structural Response Recorders to validate their use in the 
analysis. 

- Study the similarity in responses obtained from Stmctural Response Recorders and 
Strong Motion Accelerographs. 

- Develop a methodology to arrive at an estimate of the peak ground acceleration at he 
recording site from responses recorded by Stmctural Response Recorders. 
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- Develop a database of peak ground accelerations observed during Indian earthquakes 
using the available observations from SMAs and the estimated values of peak ground 
accelerations from SRRs. 

Compare the attenuation of peak ground acceleration from different events with 
respect to the various attenuation relationships available in the hterature, which are 
developed for different tectonic regions. 

Study the attenuation characteristics of peak ground accelerations recorded for 
Himalayan earthquakes for different geological regions in India and develop new 
attenuation relationships for peak round acceleration incorporating the data gmerated 
from the current study. 

Following are the salient conclusions from the current study 

- The SRR responses can be considered to come from the same sample space as that of 
Strong Motion Accelerograms. 

- The methodology of using amplification factors of resultant spectral accrieration 
from Strong Motion Accelerograms to calculate the predicted values of PGA from SA 
values of SRRs along with the usage of correlation of responses with respect to PGA 
as weighing function seems to yield good results. The predicted peak ground 
acceleration using the methodology compares very well with the recorded PGA from 
Strong Motion Accelerographs. 

- The amplification factors of the spectral acceleration for different natural periods of 
the structure indicated a decreasing trend for some natural periods and an increasing 
trend for some other natural periods for the same earthquake. The relationships for 
spectral acceleration available in the literature predict the relation of amplification 
factor with distance as one with a very weak negative dependence or a strong positive 
dependence. However, the trend predicted by these relationships is consistent for 
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periods 0.4sec, 0.75sec and 1.25 sec, which contradicts ■with the observations made 
for Indian earthquakes. Moreover, the correlation between epicentral distance and 
amplification factors was found very low. So, in the analysis of SRRs the effect of 
distance on amplification factors was not considered and it was taken as constant for 
each natural period. For SA corresponding to higher natural periods, the value of 
amplification factor with respect to epicentral distance is almost constant. 

- The average of ratios of spectral acceleration between damping ratios of 5% and 10% 
critical damping for SRRs was found to lie between 1.4 to 1.5. 

- The earthquakes in Himalayas can be classified into three different geographic 
regions with four different tectonic characteristics. The regions are Central 
Himalayas, Indo-Gangetic plains and North-East India. The earthquakes in North- 
East India can be divided into two categories. They are North East Indian Non- 
subduction zone earthquakes and North East Indian subduction zone earthquakes. 

- Most of the available attenuation relationships underestimate the accelerations firom 
Indian earthquakes. Attenuation relationships by Campbell, 1997 and Sadigh, 1997 
significantly under-predict the ground motion (in terms of PGA) for Indian 
earthquakes. Boore et al. 1997 relationship with assumed values of shear wave 
velocity gives somewhat reasonable estimates of acceleration except for Indo- 
Gangetic region andNE India. 

Youngs et al, 1997, Crouse, 1990 equations underestimate the motions from 
subduction zone earthquakes in North-East India. 

- Central Himalayan earthquakes produce the weakest motion while the Bihar-Nepal 
earthquake in the Indo-Gangetic plains has the strongest motion. 

- The rate of attenuation with distance is lowest for the subduction zone earthquakes in 
North-East India and is highest for earthquakes in Indo-Gangetic plain. 
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The relationships derived for different regions in India can be summarised as follows: 

Central Himalayan earthquakes 

ln(PG.4) = -4.768 + 0.5 86M-0:0032i?- 0.481 ln(i?); £r = 0.597 

Non-Subduction earthquakes in North-East India 

In(PG^) = -3.441 + 0.706il/- 0.828 ln(7?); c7 = 0.437 

Subduction earthquake in North-East India 

In(PG^) = -0.42 - 0.004i? - 0.241 ln(i?); cr = 0.578 

Bihar-Nepal Earthquake 

In(PG^) = 2.103 -0.006i?-0.761n(/?); cr = 0.696 
Due to rather limited data, the proposed relationships may be used only as first 
approximation for estimating ground motion. 

Based on the present work the following areas are identified with scope of research and 
refinement. 

- In the current study, only attenuation of peak ground acceleration was considered. To 
arrive at site dependant spectra, information about attenuation of spectral acceleration 
with respect to distance is very useful. The current work can be extended to derive the 
attenuation characteristics of the region in terms of spectral acceleration. The 
applicability and methodology of inclusion of recordings from Structural Response 
Recorders to the database of spectral accelerations computed from SMAs can also be 

looked into. 

- No segregation of amplification factor based on local site conditions was done. As the 
information on the local site conditions of SMA & SRR sites becomes available, an 
effort can be carried out to verify the assumption of zero correlation between 
amplification factor and site conditions. 
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The effect of epicentral distance on amplification factors of spectral acceleration was 
not considered in the present study because it showed a trend not in line with the 
reported characteristics in the literature. When data of SMA from recent earthquakes 
becomes available, these characteristics can be verified. 

Method of least squares is used to carry out the regression analysis in the current 
study. However, most of the recent attenuation relationships make use of two-stage 
regression analysis or method of maximum likelihood for carrying out the regression 
analysis. These methods are reported to be more robust. When more data from the 
Indian earthquakes becomes available, these methods can be used to develop the 
attenuation relationships. 

- A simple regression model was used in the current study to characterise the 
attenuation of acceleration. This may not give good estimate of acceleration in the 
near source region. To account for the data that will be obtained from near source 
region from future earthquakes and saturation of acceleration at near source region, a 
modified regression model can be used. 

- As the fault parameters of the earthquakes become available, the distance terms used 
can be modified to account for the fault dimensions. 


• • • 
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Appendix A 


Table Al(a): Summary of strong motion data NE-India earthquake SeptlO, 1986 


SI. 

No. 

Location 

Component 

Derived Maximum 
Peak Ground 
Acceleration In 
(cm/sec^) 

Derived Maximum 
Peak Ground 
Velocity (mm/sec) 

Derived Maximum 
Peak Ground 
Displacement (mm) 

■ 

Baithalangso 1 

omiQii 

44.51 

24.58 

41.25 

20.58 

10.91 

12.07 

3.55 

4.74 

2.73 

■ 

Dauki 

L-S 72 E 
V-VERT 

T-S 8 W 

87.61 

31.20 

88.65 

32.40 

14.14 

37.33 

3.00 

5.92 

2.56 

1 

Khliehriat 


30.30 

1.45 

45.01 

11.30 

7.71 

18.41 

2.00 

1.95 

3.86 


Mongkhlaw 


53.87 

33.55 

90.96 

33.09 

14.11 

48.08 

9.32 

6.83 

5.47 

1 

Nongpoh 


52.87 

32.86 

54.53 

21.13 

11.25 

11.32 

4.26 

2.64 

4.59 

1 

Nongstoin 

L-N 65 E 
V-VERT 

T-5 25 E 

ROi 

8.07 

13.59 

10.04 

5.95 

6.48 

4.90 

1.99 

1.66 


Panimur 

L-N 65 E 
V-VERT 

T-S 25 F 

38.35 

22.80 

47.71 

9.37 

5.61 

21.15 

2.02 

1.40 

2.65 

8 

PynuTsla 


90.95 

29.68 

74.22 

26.75 

9.28 

20.37 

324 

3.80 

4.55 

9 

Saitsama 

i3Si9BB 

110.87 

60.72 

135.86 

36.69 

20.06 

58.73 

4.90 

322 

4.92 

10 

Ummulong 



26.53 

14.25 

12.10 

238 

321 

1.56 

11 

Umrongso 

L-S 27 W 
V-VERT 

T-N 63 W 

26.73 

13.86 

31.36 

9.68 

7.81 

11.84 

130 1 

226 

4.41 

12 

Umsning 


99.49 

Al.Zl 

74.89 

28.72 

10.88 

25.67 

3.42 

3.75 

3.93 
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Table Al(b): Summary of strong motion data NE-India earthquake May 18, 1987 


SI. 

Location 

Component 

O 

Derived 

- — — 7 

Derived Maximum 

Derived Maximum 

No. 



Maximum Peak 

Peak Ground 

Peak Ground 



- 

Ground 

Acceleration In 
(ciiL*'sec^) 

Velocity (mm/sec) 

Displacement 

(mm) 
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Table A 1(c): Summary of strong motion data NE-India earthquake Feb.6, 1988 


SI. Location 
No. 


1 Baigao 


2 Baithalangso 


3 Bamungao 


Dauki 


5 Gunjung 


6 Haflong 


7 Hatikhali 


Katakhal 


9 I Khiliehriat 


10 I Mawphlang 


11 I Nongkhlaw 


12 I Nongpoh 


13 I Pynursla 


Component 


Derived Maximum 


Derived Max 


L-N 15 
V-VERT 
T-S 75 


L-N 10 
V-VERT 
T-S 80 




L-S 89 
V-VERT 
T-S 1 


* 21.43 

14.14 

3.9 

9.61 

11.27 

3.26 

23.94 

20.64 

6.9 

29.59 

13.65 

3.81 

13.7 

7.44 

1.76 

21.78 

12.26 

3.33 

16.04 

5.08 

1.3 

9.84 

6.23 

2.33 

13.05 

10.35 

2.3 

25.99 

18.01 

2.5 

30.2 

8.57 

3.41 

37.93 

18.37 

4.41 

35.7 

29.05 

6.26 

19.14 

12.87 

3.3 

36.25 

24.44 

3.87 

33.97 

19.42 

3.13 

7.88 

7.22 

1.99 

26.82 

14.02 

3.34 

23.18 

8.78 

3.18 

20.95 

7.29 

1.98 

24.89 

9.87 

3.84 

9.2 1 

10.64 

3.04 

9.65 

9.88 

4.0 

8.34 

7.75 

1.6 

78.19 

47.46 

8.45 ^ 

28.93 

15.91 

2.94 

64.86 

27.87 

3.44 


80 E 


79.61 

35.27 

57.68 


105.44 

100.76 

112.09 

26.9 

37.5 

84.62 


48.75 

15.0 

30.42 


21.24 

4.99 

9.43 1 

3.63 

11.95 

4.93 


16.06 

5.16 

14.42 

5.19 

15.39 

7.86 


8 1 Umsning 
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Table Al(d): Summary of strong motion data for Aug 6, 1988 NE-India earthquake. 
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18 

Khliehriat 

19 

Koomber 

20 

Loharghat 

21 

Mawkyrwat 

22 

Mawpblang 

23 

Mawsynram 

24 

Nongkhlaw 

25 

Nongstoin 

26 

Panimur 

27 

Pynursla 

1 

28 

Saitsama 

29 

Shillong 

30 

Silchar 

31 

Ummulong 

32 

Umrongso ’ 

33 

Umsning 


L-S 58 W 
V-VERT 
T-N 32 W 


L-N 80 E 
V-VERT 
T-S 10 E 


L-N 65 E 
V-VERT 
T-S 25 E 


L-N 65 E 
V-VERT 
T-S 25 E 


L-N 59 E 
V-VERT 
T-S 31 E 


L-N 60 E 
V-VERT 
T-S 30 E 


104.75 


83.63 

33.17 

63.3 


135.35 

81.62 

143.47 


53.23 

40.84 

50.67 


165.44 

71.62 

122.41 


207.04 

96.73 

228.28 


73.5 

35.12 

56.16 


63.18 

24.76 

89.38 


98.14 

60.8 

149.76 


76.91 

43.67 

79.04 


33.27 

69.78 

50.16 


33.67 

21.45 

28.61 

8.29 

5.89 

5.47 

64.5 

6.02 

34.17 

9.46 

55.66 

11.07 


55.19 

20.85 

44.38 


42.13 

31.55 

32.24 


87.31 

33.77 

103.63 


20.44 

9.75 

17.55 


71.52 

29.98 

98.76 


30.76 

27.33 

43.95 


51.66 

17.48 

39.98 


41.46 

38.31 

51.94 


5.29 

6.14 

6.06 


5.92 

6.35 

7.22 


7.83 

4.82 

33.89 


3.24 

3.31 

5.86 


21.1 

8.62 

21.13 


7.86 

11.28 

12.8 


6 
6 
7.76 


7.88 

2.5 

0.12 
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Table Al(e): Sui 




otion data Jan 10, 










1990 NE-India 





earthquake . 


Derived Maximum 
Peak Ground 
^ Displacement 
(mm) 


3.21 

2.02 

2.20 

1A6 

3.97 

6.00 

2.70 

3.28 

4.33 

10.75 

5.23 

16.17 

4.89 

4.57 

9.98 


3.06 

1.85 

5.87 


2.61 

2.69 

3.31 


2.94 

2.09 

6.44 

3.89 
2.48 
2,68 


6.41 
1.86 

! 4.61 


3j 
2.49 
3.11 
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Table Al(f): Summary of strong motion data for October 20,1991 
earthquake. 


SI. Location Component Derived Maximum Derived 

" Peak Ground Maximum Peak 

Acceleration In Ground Velocity 
(cm/sec^) (mm/sec) 


Uttarkashi 

























































Table Al(g): Summary of strong motion data for march 29,1999 Chamoli (U.P. hills) 


earthquake . 


Location 

Component 

Almora 

L-N53W 

V-VERT 

T-N37E 

Barkot 

L-N lOE 
V-VERT 
T-N80E 

Chinyalisaur 

L-N43E 

V-VERT 

T-N47W 

Ghansiali 

L-NOOE 

V-VERT 

T-N90E 

Gopeshwar 

L-N 70 W 

V-VERT 

T-N20E 

Joshimath 

L-N80E 

V-VERT 

T-NIOW 

Lansdown 

L-N70E 

V-VERT 

T-N20W 

Roorke 

L-N55W 

V-VERT 

T-N35E 

Tehri 

L-N63W 

V-VERT 

T-N27E 

Ukhimath 

L-N 15E 
V-VERT 
T-N75W 

Uttarkashi 

L-N 15 W 

V-VERT 

T-N75E 


Denved Maximum 
Peak Ground 
Acceleration In 
(cm/sec“) 

Derived 

Maximum Peak 
Ground Velocity 
(mm/sec) 

Derived 
Maximum Peak 
Ground 
Displacement 
(mm) 

26.17 

23.50 

3.70 

26.89 

12.70 

0.90 

27.74 

18.90 

3.00 

16.88 

8.30 

0.40 

19.59 

7.70 

0.30 

22.26 

12.40 

0.50 

50.65 

31.20 

2.50 

48.64 

31.40 

2.70 

43.92 

33.10 

3.10 

71.42 

33.10 

5.00 

38.74 

20.00 

5.30 

81.84 

40.80 

0.70 

195.08 

225.50 

52.20 

153.73 

75.00 

20.90 

352.83 

453.10 

122.80 


69.55 

40.45 

6228 

5.24 

11.01 

6.32 

5539 

17.25 

45.89 


53.43 

33.54 

61.06 


89.06 

46.45 

94.71 


52.74 

22.22 

62.32 


31.80 

27.20 
89.40 

1.80 

2.50 

2.10 

42.80 

12.20 

32.90 


46.0' 

21.1 

53.8 


68.50 

40.50 

54.50 


35.50 

15.80 

46-70 


SS 
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Table Al(h): Summary of strong motion data, March 29,1999 Kangra (U.P. 
earthquake. 
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Table A2 (b): SMA stations installed in Uttarpradesh array and their geographical 

coordinates 


SI No: 


Station 


Longitude 


Latitude 


Almora 

Ba geswar 

Barkot 

Bhatwari 

Chak 

Chamoli 

Chau 

Chinsyali 

Didihat 


10 Gwal 

1 1 Haldwani 

12 Kapk 


13 Kamapra yag 

14 Kosani 


15 Koteswar 


Koti 

Lansdowne 


18 Nainital 


19 Pithoearh 


20 Purola 


21 Ramnagar 
~72 Ranikhet~ 
Rati 


Roorkee 


25 1 Srinagar 


26 Uttarklashi 



Table A2(c): SMA stations installed in Kangra arrc^ and their geographical 
coordinates 
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Dalhousie 


D h armpur 
Dh armsala 
Drang , 
Chumarwin 


Jawali 


Jhatingri 


Jhungi 


Kan 


Kotla 


Kulu 

Kunihar 


Ladbharol 


Lambagaon 


Nadaun 


33 Nagrota Bagwan 

34 Namhol 


35 


36 Pandoh 


37 Patlander 

38 Pragpur 


Rakh 


Rewalsar 


41 I Rirkmar 


42 Sandhu 


43 


44 Shahpur 


45 Shimla 


46 Sihunta 


47 Sundemagar 


48 Siindla 


49 Sunni 


50 I Surani 


76 

11 

77 

22 

76 

44 

76 

11 

77 

13 

76 

5 

76 

54 

76 

2 





























































































































Appendix B 

Table Bl(a): Recorded SRR values and the predicted PGA ’sfor August 6, 1988 N-E 
India earthquake 
























































































































































































































































































































































































































































































Appendix B 


Table Bl(b): Recorded 

Bihar-Nei 

I SL 1 Stations I 


SRR values and the predicted PGA ’s for August 

jal earthquake. 

dist. Accelerations in 

5% damping ” 10% damping T” 

0.4 I 0.75 I 1-25 ^4 I 0.75 I 1.25 




21 , 1988 , 



























































































































































































































































Appendix B 
Table Bl(p 


Bl(c): Recorded SRR values and the predicted PGA ’s for October 20, 1990, 
Uttarkashi arthquake. 


Stations 


accelerations % 


Au t 

Barnath_ 
Bajaiira 
3arjar_ _ __ 

Biiaspur 


Dharamashala 


Drane^ 


Ghuinanvin 


ur 


Jottindrana^ar 


Jwalamukhi 


era 


Kasull 


Kulu 


Manall 


Naduan 


Naban 


Narkanda 


Pandoh 


Parwaiioo 


ooli 


0.4 


_15 0. 0 2.83 
2^4.2 6.67 


0.4 0.75 


2.49 1 0.99 


1.28 


8.06 i 6.01 


Rohm 


Shimla 


Sirka 




Una 

Ambala 


Faridabad. 


Gohana 


Jagadhari 

JhaJJar 

Jind 


Kaittial 


Ktirukst,ietra 


Narwana 


Rohtak 


Sohna 


8.21 


3.27 


241.8 3.58 


1.65 


1.03 


1.18 0.76 


2.39 I 1.72 I 8. 


2.34 2. 


1.23 2.76 


2.71 2.33 3.08 


1.59 I 1.12 1 6.16 


159.71 4.85 I 2.14 1 2.01 


285.71 2.18 


171.2 


211.3 


1.94 2.21 


116.8 11.82 14.09 


255.6 8.43 2.27 


197.1 12.54 3.25 


197.ll 6.18 1 3.06 


2.69 


0.78 


4.26 


197.1 8. 


320.6 2. 


261.4 3.32 


189.9 3.30 1.18 ! 0.34 


178.6 9.01 1.00 


3.82 1.76 


8.14 2.64 


4.76 


112.61 15.71 I 6.25 I 2.61 


218.ll 1.84 


0.59 


1.86 


1.04 2.13 


3.19 7.76 


2.69 4.99 


0.45 1.06 


8.10 7.66 


6.86 1.62 


1.15 4.65 2.94 


1.55 I 5.53 I 2.65 


2.23 


1.09 


1.30 I 2.91 


2.91 


244.9 2.91 


199.8 5.49 


249.7 3.83 

194.2 7.48 


287.5 3.64 


299.4 3.08 


289.8 4.91 


1.54 


0.73 2.05 


n I 4.11 
93 


1.91 


Mwsm 


11 I 0.73 

sTl \J2 


1.01 


1.86 
1.85 I 2.91 


269.5 5.90 2.37 3.25 4.21 

154.61 10.47 I 5.76 \.12 j 6.02 

1.67 0.58 

1 2.00 1.45 1 2.57 


310.8 1.98 1.67 0.58 

279.3 2.91 2.00 1.45 

244.9' 3.82 2.03 0.90 

227.3 4.50 6.29 2.82 

284.9 5.42 2.23 0.78 

228.1 5.98 1.38 1.13 

299.3 3.18 0.52 0.53 


1.25 


0.85 


0.77 


1.66 

1.35 



1.31 


0.98 


1.51 


2.05 


1.16 


1.84 0.94 


321.91 2. 


a 


1.92 2.01 


1.18 

L70 

iff ~ 

3.71 ■■ 

3.9 r~ 

1.19 


1.16 


2.86 


1.76 


1.16 


1.90 

1.93 


2.12 


0.83 


2.48 


1.16 


3.21 


2.97 


0.92 


5.81 


2.60 


2.81 


2.57 


3.25 


0.88 


2.23 


0.88 


2.21 


1.60 


2.67 


4.12 


5.03 


1.18 

1.09 


1.55 


1.32 


2.42 


1.45 


1.54 


1.83 


2.98 


2.98 


0.97 


1.38 


1.63 


3.12 


1.59 


1.49 


0.81 


1.48 
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Afa?a!earh 

A lmor a 

Am roha 

As k ot 

Ba daun 

Badrinath 


59 


60 


6 1 1 Bhatw’ari 


62 Rijnor 

63 Bilaspur 

64 Champawat 
Chandausi 

66 Chandpur 

67 Dakpathar 

Dehr aduo 

69 [ Deoband 


70 


71 Dharasu 


72 Dharchula 


73 Didihat 


74 


75 I Haldwani 


76 


77 Hardwar 


78 Joshimath 


79 Kaima 


80 


81 


82 


83 Kausani 


84 Koti 


85 Lansdowne 


86 Meenit 


87 Moradabad 

88 Mussorie 


1S4.8 


299,2 1 2.76 


6.16 


5.35 


1.24 .1.62 


1.42 


4.32 


2.18 


4.97 


84.8 


152.2 


164.6 


51.7 

d 

156.7 


177.3 





a 


Nainaba 


Naiinital 


Narendra 


Paurt 


Pilibhit 


99 


100 1 Ranikhel 


101 Rishikesh 

102 Roorkee 


103 


104 

105 1 Sambhal 


7.61 


20.56 


2.05 


7.12 


4.15 


3.25 


2.63 


7.32 


3.48 


3.49 


4.26 


2.33 


3.13 


2.49 


20.21 


16.00 


4.11 


3.04 


16.81 


4.11 


8.73 


24.43 


6.98 


5.97 


1.49 


4.89 


17.88 


1.99 


5.98 


8.21 


2.62 


1.25 


4.36 3.10 

3.88 1.41 

4.92 5.31 

1.45 1.19 

17.45 13.96 


3.79 I 1.45 


3.13 1.65 


2.18 0.74 


0.82 


0.92 


1.90 1.75 


2.69 7.92 


1.98 0.70 


2.83 0.78 


2.36 0.93 


11.64 4.46 


5.88 3.99 


0.79 2.36 


2.49 1.70 


2.14 0.77 


2.80 


2.72 


3.64 


0.94 


0-58 


11.25 


1.10 


0.73 


0.98 


3.42 


1.34 


1.13 


0.45 

6.82 


1.18 


1.28 


1.10 


9.96 


2.14 


2.33 


11.45 


2.99 


3.06 


0.97 


1.77 


5.03 


1.35 


2.09 


2.09 


1.50 


5.81 


3.61 


3.09 1.10 

2.40 I 6.23 5.53 









































































































































































































































































































































Table Bl(d): Recorded SRR values and the predicted PGA^s for March, 29 1999, 
Chamoli earthquake. 


SI. 

No. 

Stations 

Dist. 

Acceleration in %g | 

(km) 

5% 1 

10% 1 

PGA 


0.4 

sec 

0.75 

sec 

1.25 

sec 

0.4 

sec 

0.75 

sec 

1.25 

sec 

1 

Barkot 

131.3 

10.19 

3.26 

1.22 

7.03 

EHII 

1.12 

2.81 

2 

Bhatwari 

81.5 

7.34 

2.55 

1.12 


■IHI 

0.71 

1.76 

3 

Dehradun 

131.3 

21.51 

BHI 

3.67 

18.65 

8.87 

3.16 

7.52 

4 

Dharasu 

118.4 

3.98 

WEm 

1.53 

1.22 

1.22 

0.92 

1.35 

5 

Gopeshwar 

4.8 

64.53 

45.16 


64.12 


20.59 

33.53 

6 

Joshimath 

26.9 

14.07 


9.99 



6.63 

8.12 

nr 

Kamprava:a 

20.6 

6.93 



5.91 


1,.T. 

3.57 

"s~' 

Kausani 

BiMM 

9.17 



8.46 


1.43 

3.18 

s 

Lansdown 

i!i!3Wi 

3.98 

EH 

0.92 

1.83 

1.02 

0.61 

1.23 

m 

Mussourie 

127.8 

19.88 

msm 

2.55 

17.13 

7.75 

2.04 

6.44 

ni 


43.6 

17.43 

IM 

2.96 



2.45 

5.22 

12 

Pauri 

108.6 


|||g|| 

4.49 

\WEEM 

■m 

4.18 

3.63 

13 

Roorkee 

158.7 

itiW 


2.85 

Mifrai 

■H 

2.14 

3.82 

14 

Rudraprayag 

44.7 

4.69 

2.04 

3.98 

4.28 

1.53 

2.75 

2.83 

15 

Tehri 

88.0 

25.08 

9.28 

3.16 

16.11 

8.56 

ins 

7.17 

16 

Uttarkashi 

111.8 

28.95 

5.81 

5.50 


5.40 

2.96 

6.77 


Table Bl(e): Recorded SRR values and the predicted PGA ’s for April 26 1986, Kangra 


earthquake. 
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Table B2: SRR stations installed and their geographical coordinates 


Station 


-Along 

j^mdila 

Chou kham 

Darin g 

Dirang 

Gcnsi 


Hanoi i 

_Itanagar _ 

10 Khonsa 

11 Oyan 

12 _Pangin 

13 ^j'iisighat 

14 Raga 


15 


16 1 Tawan 


17 Yemboun 


18 


19 


20 


21 

22 


23 


24 Bihpuria 


25 


26 Bijoynagar 

^ 27 Bilasip ara 

28 BokaJkrat 


29 


30 I Chaparmukh 


31 


32 


33 


34 


35 I Di-iin 


36 I Dibruearh 


37 

38 


39 I Doomdooma 


40 Dudhnoi 


41 I Dullabchera 

42 


43 


44 


45 


46 Goriipur 

47 Guwahati 


48 Haflong 

49 Hailakandi 


50 Hojai 


Longitude 


Latitude 


Deg. 1 Min. 


Aru nac hal 

Arunachal 


Arunachal 


Arunachal 

Arunachal 


Arunac 


Arunachal 

Aiunach:^ 

Ar una chal 

Arunachal 

Arunachal 

Arun achal 

Aru na chal 

Arunachal 

Arunachal 

Arunachal 

Arunachal 

Assam 


Assam 


Assam 


Assam 


Assam 


Assam 


Assam 


Assam 


Assam 


Assam 


Assam 


Assam 


Assam 


Assam 


Assam 


Assam 


Assam 


Assam 


52 

27 

35 

12 

27 

28 

38 

26 

21 

31 

26 

49 
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Station 


Jamugiri 

^ Jorhat __ ^ 
Kamargaqn_^ ____ 

Kamp ur 

, _Karimganj 

j K ollabar 

58 Kukichera 

59 La khi pur 

60 Lalagiiat^ 

Lank a 

^ 62 ^ ' "“redo 

63 Likh apani 

64 Lumding 

65 _Nlairabari 

66 M a kum 

67_ Mag aldai 

68 Mariani 


69 Moranhat 


70 Naharkatiya 


Nalbari 


Nazira 


Neamati 


74 1 Nilbaean 


75 Nogora 

76 North-l akhimpur 


77 Paneri 


78 


79 Rangia 

80 Sarupeta 


81 


82 


Silchar 


85 

86 I Subankrata 

87 


88 Tinsukia 


89 I Titabar 

90 


91 I Ararira 


92 

93 Bairagania 

94 Barauni 


95 I Bettiah 


Assam 
A ssam 
Assam _ 

Assa m 

Assam 

A ssam 

Assam 

Assam 

Assam 

Ass am 

Assam 

Assam 

Assam 

Assa m 

Assam 

Assam 

Assam _ 

Assam 


Assam 


Assam 


Assam 


Assam 


Assam 


Assam 


Assam 


Assam 


Assam 


Assam 


Assam 


Assam 


Assam 


Assam 


Assam 


Assam 


Assam 


Assam 


Assam 


Assam 


Assam 


Assam 


Bihar 


Bihar 

Bihar 


Bihar 


Bihar 


Bihar 


Bihar 


Bihaj;_ 

"Bihar" 


Bihar 


27 


25 


26 


27 



26 


26 


26 


25 


26 


27 


26 


26 


26 


26 


27 


26 


24 


26 
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Station 


Latituie 


Jwalamuki 


Kandaghat 


1531 Kan 


Kasuali 


156 1 Khoksar 


1571 Kulu 


Manali 


1601 Mandi 


161 Nachar 


162 Naduan 


Nahan 


1651 Narkanda 


166 


167 Palawur 

168 Pandoh 


1691 Parwanoo 


Pooh 


Rohm 


Sirkhaghat 


174 


1751 Sundamagar 


176 Una 


177 Talwara 


178 Ambala 


179 Bahaurgarh 


180l Ballabhgarh 


181 1 Faridabad 


182 1 Gohana 


183 


184 I Jagadhari 


1851 Jhajjar 


186 Jind 


187 Kaithal 


1881 Kumksihetra 


1891 Narwana 


Himachal 


Himachal 


Himachal 


Himachal 


Himachal 


Hiniachal 


Himachal 


Himachal 


Himachal 


Himachal 


Himachal 


Himachal 


Himachal 


Himachal 


Himachal 


Himachal 


Himachal 


Himachal 


Himachal 


Himachal 


Himachal 


Himachal 


Himachal 


Himachal 


Himachal 


Himachal 


Pimiab 




Rohtak 


Sohna 


1951 Cherapunji 


196 Dalu 


197 Dauki 


198 Jowai 


199 1 Mawplang 
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SI. Station 

No. 


25 1 Belonia 

25,2 Bishalgarh 

253 Kailashahar 

254 Kamalpur 

255 Kiiowai 

256 Nutarbari 

257 Radhakishorepur 

258 Sabrum 

259 Afazalgarh 

260 Almora 

261 Amroha 

262 Askot 

263 Badaun 


2641 Badnnath 


265 


2661 Baharaich 


267 Balrampur 

268 Banbassa 


269 


Latitude 


Deg. I Min. 


15 



Bhatwari 


274 


275 Bulandshahar 


276 Chawawat 


2771 Chandausi 


2781 Chandpur 


2791 Dakpathar 


2801 Dehradun 


281 1 Deoband 


2821 Dhamapur 


Dharasu 


Dharchula 


Didihat 


Domariyaganj 

Gonda 


2881 Gopeshwar 


2891 Haldwani 


290 


Hardwar 


Jarwal-road 


Joshimath 

2941 Kairana 

295 Kaisarganj 

296 Karanprayag 

297 Kashipur 


298 Kausani 

299 Lansdowne 


3001 Mthura 


29 

10 

27 

15 

30 

16 

28 

12 

29 

30 

29 

48 

27 

28 
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Station 


Latitude 


301 Meerut 

302 Moradabad 

303 Mussorie 

304 Muzzafama^ar 

305 Nagina 

306 Nainabag 

307 Nainital 

308 Nanpara 

309 Narora 

310 Nautanwa 

311 Nepalganjroad 

312 Nicblaul 

313 Pauri 

314 Pilibhit 

315 Pithoragarh 
316 1 Puranpur 

317 

318 Raniketh 

319 Rishikesh 


320 Saharanpur 


Sambhal 


Shamli 


3231 Teiuu 


324 


325 Thai 


326 Uska-bazar 


327 Uttarkasi 


328 Alipurdwar 


329 1 Belakoba 



Haldibari 


Kurseon 


Madarihat 


336 

337 1 Ramshai 

338 Siliffliri 

339 Champa 

340 Pathankot 

341 Sh imla 

342 

343 1 Forbesganj 
344 Marghetia 


Delhi 

Gorakhpur 

Kalagarh 


3481 Shillon 


Vxharali 

Koti 


32 

16 

31 

8 

31 

49 

26 

18 

27 1 

17 

28 

40 

26 

45 

29 

30 

25 

. 34 

26 

45 ^ 

30 

35 






























































































































